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in ■ ABSTRACT 

o 
o 

■ The effect of smooth cloud boundaries on the interaction of steady planar shock 
^ ' waves with interstellar clouds is studied using a high-resolution local adaptive mesh re- 

. finement technique with a second-order accurate axisymmetric Godunov hydrodynamic 

scheme. A three-dimensional calculation is also done to confirm the results of the two- 
dimensional ones. We consider an initially spherical cloud whose density distribution is 
. ^ flat near the cloud center and has a power-law profile in the cloud envelope. Our model 

\^ can be specified by three model parameters: the power law index of the density profile 

(n), the density contrast (x), and the Mach number of an incident shock relative to the 
sound speed of the preshock ambient medium (M). The evolution of the shocked cloud 
is divided into four stages: (1) the initial transient stage, in which the initial shock 
. is transmitted to the cloud, and at the same time is refiected at the cloud surface to 

form a bow shock or bow wave; (2) the shock compression stage, in which the shocked 
Qh| cloud is compressed in the direction of shock propagation; (3) the postshock expansion 

O I stage, in which the cloud expands in the radial direction; and (4) the cloud destruction 

' stage, in which the shocked cloud is shredded into small fragments and is eventually 

■ destroyed. This evolution is qualitatively consistent with previous simulations of shocks 
interacting with uniform clouds with sharp cloud boundaries. However, smooth cloud 

^ . boundaries significantly affect the morphology of the shocked clouds and cloud destruc- 

tion timescales. When an incident shock is transmitted into a smooth cloud, velocity 
gradients in the cloud envelope steepen the smooth density profile at the upstream side, 
resulting in a sharp density jump having an arc-like shape. We refer to this density 
jump as a "slip surface" because the velocity is sheared parallel to its surface. Such a 
slip surface forms immediately when a shock strikes a cloud with a sharp boundary. For 
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smoother boundaries, the formation of shp surface and therefore the onset of hydrody- 
namic instabihties are delayed. Since the slip surface is subject to the Kelvin-Hclmholtz 
and Rayleigh- Taylor instabilities, the shocked cloud is eventually destroyed in 3 — 10 
cloud crushing times, where the cloud crushing time is the characteristic time for the 
shock to pass through the cloud. The ratio of the cloud destruction time to the cloud 
crushing time depends mainly on the smoothness of the density profile and only weakly 
on the density contrast for strong shocks. We construct analytic models of cloud drag 
and vorticity generation that compare well with the numerical results. After complete 
cloud destruction (in several cloud crushing times), small blobs formed by fragmenta- 
tion due to hydrodynamic instabilities have significant velocity dispersions of the order 
of 0.1 vj,, where Vb is the shock velocity in the ambient medium. Our simulations show 
that, irrespective of the smoothness of the initial density profile, it is difficult to sustain 
the vortical motions inside the cloud because the cloud is being destroyed by hydro- 
dynamic instabilities. This suggests that turbulent motions generated by shock-cloud 
interaction arc directly associated with cloud destruction. The interaction of a shock 
with a cold H I cloud should lead to the production of a spray of small H I shreds, which 
could be related to the small cold clouds recently observed by Stanimirovic &; Heiles 
(2005). We also calculate the dependence of the velocity dispersion (Av) on region size 
(R), the so-called linewidth-size relation, based on a three-dimensional simulation. The 
lincwidth-size relation obtained from our simulation is found to be time-dependent. In 
the early stages of cloud destruction, the small-scale fluctuations dominate because of 
the nonlinear growth of the Kelvin-Helmholtz instability, and thus the linewidth-size 
relation is more or less flat. In the later stages, the small-scale fluctuations tend to 
damp, leading to a linewidth that increases with size. A possibility for gravitational 
instability triggered by shock compression is also discussed. Only radiative shocks with 
7 < 4/3 can lead to induced gravitational collapse. We show that in the absence of 
significant nonthermal motions behind the shock, star formation can be induced in an 
initially uniform cloud by a radiative shock only if it is not too strong. Our results in- 
dicate that the postshock gas has significant nonthermal motions, and this can further 
limit the onset of gravitational collapse. 

Subject headings: hydrodynamics — ISM: clouds — shock waves — supernova remnants 
— turbulence 



1. Introduction 

The interaction of shock waves with interstellar clouds is a very common event in interstellar 
medium (ISM). For example, supernova explosions produce large blast waves that propagate into 
the ISM. The blast waves strike interstellar clouds, compressing them and destroying them. This 
process plays an important role in the formation of multi-phase structures in the ISM (Cox & 
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Smith 1974; McKee &; Ostriker 1977). Recently, several authors have proposed it as a mechanism 
for generating the supersonic turbulence observed in the ISM (e.g., Kornrcich & Scalo 2000; MacLow 
& Klessen 2003; Vazquez-Semadeni et al. 2000). Compression by shocks may be able to trigger 
star formation in interstellar clouds (Elmegreen & Lada 1977). In the early universe, supernova 
explosions of the first stars may be responsible for the first significant generation of turbulence 
and the first metal enrichment in the intergalactic medium after the Big Bang. Other energetic 
events, such as stellar winds of massive stars, bipolar oTitflows of young stellar objects, cloud- 
cloud collisions, and spiral density waves can also produce interstellar shocks that interact with the 
surrounding clumpy ISM. Hence, understanding how shock waves interact with interstellar clouds 
is an important step in understanding the structure and evolution of the ISM. 

Many authors have studied the interaction of shock waves with interstellar clouds, either 
analytically (e.g., McKee & Cowie 1975; Spitzer 1982; Heathcote k Brand 1983; McKee et al. 1987; 
Nulsen 1982) or numerically (Sgro 1975; Woodward 1976; Nittman et al. 1982; Tenorio-Tagle 
& Rozyczka 1986; Bedogni & Woodward 1990; Stone & Norman 1992; Klein, McKee & Colella 
1990, 1994 (hereafter KMC), Mac Low et al. 1994; Dai k Woodward 1995; Xu k Stone 1995). 
Such calculations are closely related to the interaction of a cloud with a wind (KMC; Schiano, 
Christiansen, k Knerr 1995; Gregori et al. 2000). Analytic studies have focused mainly on the 
structure and evolution of a shock transmitted into a cloud before cloud destruction. After the 
shock passes through the cloud, hydrodynamic instabilities such as Rayleigh- Taylor, Richtmyer- 
Meshkov, and Kelvin-Helmholtz begin to develop and may eventually destroy the entire cloud. 
It is difficult to follow such highly nonlinear stages of evolution analytically, and thus numerical 
calculations are needed. KMC performed comprehensive two-dimensional numerical studies of the 
interaction of shocks with uniform clouds based on the adaptive mesh refinement (AMR) technique. 
They showed that small, nonradiative clouds are destroyed in several cloud crushing times, where 
the cloud crushing time is the timescale for a shock to propagate through the cloud (see §2 for 
the definition). They also demonstrated that good spatial resolution (at least 10^ cells per cloud 
radius) is needed in numerical calculations in order to follow the cloud destruction accurately. Stone 
k Norman (1992), Klein k McKee (1994), and Xu k Stone (1995) performed three-dimensional 
numerical simulations on this problem, showing that overall evolution of the cloud destruction is in 
good agreement with the two-dimensional results of KMC, although vortex rings observed in the 
two-dimensional calculations are unstable in three dimensions and the cloud destruction proceeds 
more violently in all directions. 

Recently, Klein et al. (2000, 2003) performed a series of Nova laser experiments investigating 
the evolution of a high density sphere embedded in a low density medium after the passage of a 
strong shock (with a Mach number of ~10). The results of these experiments can be compared 
directly with interstellar shock-cloud interaction under certain scale transformations (Ryutov et 
al. 1999). The comparison of their experimental results with two-dimensional simulations can 
reproduce the morphological evolution of the experimental shocked spheres reasonably well. In later 
stages of cloud destruction, the three-dimensional instability of vortex rings plays a significant role. 
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and fully three-dimensional calculations are required to follow the evolution. Follow-up experiments 
on the Omega laser confirmed the presence of vortex ring instabilities, which were shown to be in 
good agreement with the predictions made by the Widnall instability (Robey et al. 2002). 

Previous numerical studies have assumed that the initial clouds have a uniform density and 
sharp boundaries. However, real interstellar clouds have internal density gradients. For example, 
in molecular cloud cores with no embedded young stars, the density distributions are flat near the 
center and decrease with the distance from the cloud center (Bacmann et al. 2000) . Observations of 
NGC 7293 (the Helix nebula) imply the presence of large radial density gradients in the cometary 
cloud cores there (Burkert & O'Dell 1998; O'Dell et al. 2000). The actual structure of the clouds 
is likely to be complex, but observations suggest that it would be worthwhile to study the effect of 
smooth density gradients on the cloud-shock interaction. Kornreich & Scalo (2000) have carried out 
the first such study. They derived an approximation for the distribution of vorticity generated in 
a shocked cloud with a smooth density profile. They argued that the vortex motions generated by 
the shock could be converted into turbulent motions observed in interstellar clouds provided that 
the motions are confined in the parent cloud. However, from their work, it was unclear whether 
such a cloud can survive after the shock passage. In fact, many numerical calculations have shown 
that uniform clouds, which have sharp density jumps at the cloud-intcrcloud boundaries, are easily 
destroyed primarily under the influence of Kelvin-Helmholtz and Rayleigh- Taylor instabilities (e.g., 
KMC). Recently, Poludnenko, Frank, &; Blackman (2002) observed complete destruction for a cloud 
with a smooth but steep density gradient. It therefore remains an open question whether shocks 
can maintain turbulent motions in interstellar clouds. 

In this paper, we study the effects of density gradients on the cloud destruction process, using 
our higher-order Godunov AMR hydrodynamic code (Truelove et al. 1998; Klein 1999). We first 
present two-dimensional axisymmetric calculations. A three-dimensional simulation is also done 
to clarify the property of turbulent motions generated by the shock-cloud interaction. In §2 we 
summarize the formulation of our problem. We show that at least 10^ grids per cloud radius are 
needed to follow the evolution of smooth clouds accurately, just as in the uniform case. An overall 
description of shocked clouds is presented in §4, in which we demonstrate that the initial smooth 
density profile steepens after shock passage, and a density discontinuity is formed at the upstream 
side of the shocked cloud. A significant velocity shear develops across the discontinuity. The velocity 
shear makes this "slip surface" unstable to the Kelvin-Helmholtz instability, leading to complete 
cloud destruction even for smooth clouds. We also investigate the dependence of cloud evolution on 
the initial model parameters such as the density profile, shock strength, equation of state, and cloud 
shape. In §5, we estimate how the cloud is accelerated by the shock, which affects the growth rate 
of the Kelvin-Helmholtz instability, and therefore the timescale of cloud destruction. In §6 and §7, 
the formation of the slip surface and the subsequent cloud destruction are examined in more detail. 
Since the entire cloud is evcntvially destroyed by the shock, the vorticity generated by the shock is 
not confined to the original cloud. Instead, the shear motions are converted into random motions 
of small fragments after cloud fragmentation. In §8, the vorticity generated by shock passage is 
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compared with a simple analytic model. In §9, we summarize our results, and discuss how density 

gradients affect the observed morphologies of shocked clouds. We discuss several properties of 
turbulent motions generated by shock-cloud interaction, computing synthesized spectra and the 
velocity-size relation. Finally, we derive the condition under which shock waves are capable of 
triggering gravitational contraction of shocked clouds resulting in star formation. 



2. Formulation of the Problem 

2.1. A Cloud with a Smooth Boundary 

In this paper, we investigate the hydrodynamic interaction of a planar shock with an interstellar 
cloud, focusing on the effect of a smooth cloud boundary. We consider a steady planar shock 
moving through the intercloud medium at a velocity Vb at infinity, corresponding to a Mach number 
M = Vb/Cio, where Qo = iriiPiO / PioY ^"^ is the adiabatic sound speed in the intercloud medium. 
The subscript i denotes the values in the intercloud medium. The subscript denotes the values 
in the preshock state. (In the following, the subscript 1 is used for the quantities in the post- 
shock state.) This shock hits a spherical cloud whose density distribution is given in a cylindrical 
coordinate system (tu, z) centered on the cloud by 

where r = (tu^ -|- z'^)^^'^, r^o is a characteristic radius, pcQ is the density at the cloud center, and piQ 
is the density at infinity. The initial density contrast between the cloud center and the intercloud 
medium is 

X = Pco/Pio ■ (2) 

This density distribution is essentially the same as Kornreich & Scalo (2000) , although they applied 
this profile to a cylindrical cloud with a two-dimensional Cartesian geometry. For comparison, we 
shall follow a couple of runs with Cartesian geometry. For n ^ oo, this model cloud reduces to 
a uniform cloud with a sharp boundary. Because we do not include the effect of gravity in our 
simulations, the cloud is assumed to be initially in pressure equilibrium with the ambient medium 
with a constant pressure Pq everywhere in front of the shock. We also focus on nonradiative shocks, 
and therefore the ratio of specific heats is set to 7 = 5/3 for most runs, although we calculate a 
couple of models with 7 = 1.1 to examine the effects of radiative cooling. 

In order to minimize the effects of the density variation along the initial shock front, we must 
begin our calculation with the shock in a region where the density variation is less than a few 
percent of the density at infinity {pio). Thus, the distance from the cloud center to the initial shock 
front is larger for clouds with smoother density gradients and/or larger density contrasts. 

To measure the degree of mass concentration in the model cloud, we define two characteristic 
radii: a core radius, Tco, and a cloud radius, rd. The core radius denotes the size of the central 
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relatively flat region and is equal to the characteristic radius, Vco- The density at Tco is {x~^ + 
l)pco/2, which approaches Pco/2 for large %. The cloud radius, rd, denotes the overall size of the 
cloud, and is defined to be the radius at which the density is twice that of the ambient medium 
(2pio), 

rci = ix- 2)^/V,o. (3) 

We assume that rd > rco, which is true for % > 3. These two characteristic radii divide the model 
cloud into two parts: the core (r < rco) and the cloud envelope (rco < r < rd). The region r > rd 
is the intercloud medium. A third characteristic radius, the rms radius, is introduced in §2.3 below. 

We define the core mass and the cloud mass in terms of density: The core mass, mco, is the 
mass of gas with p > 0.5(1 +x~-'^)pcO) so that initially it is the mass inside rco- The cloud mass is the 
mass with p > 2pii, so that initially it is the mass inside rd- These masses begin with values mco,o 
and md^O) respectively, and then decrease due to shock heating and to mixing with the shocked 
intercloud medium. The degree of mass concentration in the initial cloud is measured by the ratio 

n^^^, (4) 

which is unity for a uniform cloud, where mdfi is the initial cloud mass. In the following, the initial 
core mass, ?tIco,Oi is defined as the mass contained within r < rco- For example, for n = 2 and large 
X, we have TZ 4x^/^ - 3vr/4. 

The thickness of the cloud boundary (/ibd) can be approximated as the density scale height at 
r = rco, 

hhd = i - — 7 ) -Tco ■ (5) 

In real interstellar clouds, the value of n is likely to vary widely. For example, in dense 
molecular cloud cores with no embedded young stars, the density profiles are flat near the center, 
follow a power-law profile of p oc in the intermediate radii, and fall steeply as p oc r~^ or even 

near the edges of the cores at r ~ 10^ AU (Bacmann et al. 2000). Observations of the Helix 
nebula imply that the density distributions in the cloud cores can be fit with power-law profiles 
with n ~ 3 - 9 (Burkert k O'Deh 1998; O'DeU et al. 2000). In the following, we vary n in the 
range of 2 < n < oo. 

We do not take into account the effect of thermal conduction, which can play a role in evolution 
of non-magnetized clouds. According to KMC, non-magnetized clouds with x ^ 10^ are ablated 
by evaporation in a time comparable to the cloud crushing time if the shocked intercloud medium 
is hot enough that the mean-free path is comparable to the cloud radius, so that the conduction is 
saturated (Cowie & McKee 1977). To neglect the effect of thermal conduction, we assume either 
that the intercloud medium is not that hot, or that there is a magnetic field that is too weak to be 
dynamically significant, but is strong enough to inhibit thermal conduction. 
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2.2. Time Scales 

Following KMC, we define two characteristic time scales governing the evolution of the shocked 
cloud: (1) the shock-crossing time, 

*sc = — ^, (6) 

Vb 

which is the time scale for the shock in the intercloud medium to sweep across the cloud, and (2) 
the cloud-crushing time, 

~ Vb MC\o MCco ^ ' 

= 0.98 X lOViMfoiC-oJs yr (8) 

where ri = rco/lpc, Miq = M/10, Cco,5 = Cco/W^ cms~^, Cco is the sound speed in the preshock 
cloud, Cio is the sound speed in the preshock ISM, and tec is the time scale for the cloud to be 
crushed by the shock transmitted into the cloud. Here, the cloud-shock velocity is approximated as 
Vg — Vb/x^^'^^ since the pressure behind the cloud shock is comparable to that behind the intercloud 
shock for a strong shock {pcoVg ~ Piovf, where Vg is the velocity of the cloud shock). We note that 
for smooth clouds, the actual mean velocity of the cloud shock is estimated as Vg — Vb/ {x)^^^^ 
where the mean density contrast, (x), is generally smaller than %. Therefore, for a smooth cloud, 
the cloud shock can reach the initial cloud center by the time of 1 tec ix) /x (< 1 ^cc)- 

According to McKee et al. (1987), the typical pressure variation timescale for a dense cloud in 
a Scdov- Taylor blast wave is given by tp ~ O.lRc/vij, where Rc is the distance from the origin of 
a blast wave. Therefore, our assumption of a steady shock requires that tec *C tp, or equivalently, 
that the cloud radius be smaller than 0.1 Rc/x^^'^ ■ This is known as the small cloud approximation 
(see §2 of KMC for more detail). 

To investigate the overall evolution of the shocked clouds, we also define the following four 
characteristic timescales: (1) the drag time, tdragi which is the time when the cloud velocity relative 
to the velocity of the post shock ambient flow has decreased by a factor of 1/e (see §5), (2) the 
maximum expansion time, tm, which is the time when the effective cloud radius in the tu-direction 
(a) has increased to 90 % of its maximum value (see §5), (3) the cloud destruction time, tdesti which 
is the time when the mass of the main axial fragment has become smaller than the initial cloud 
mass by a factor of 1/e (see §7), and (4) the mixing time, tmix, which is the time when the mass of 
the region with p > 2 pn has reached half the initial cloud mass (see §7). 



2.3. Global Quantities 

Following KMC, we define several integrated physical quantities to study the cloud evolution. 
In KMC, a two-fiuid treatment that followed both the cloud and intercloud matter allowed the 
global quantities to be evaluated by integrating the physical quantities only in the cloud material. 



-8- 



In this paper, we use a density threshold, pth; to determine which material can be considered to be 
"cloud material." The corresponding threshold mass is 



mth 



I pdV . 



(9) 



The threshold density /9th is chosen such that mth corresponds to either the characteristic initial 
core mass (mco,o = 1^ '''co, t = 0]) or the initial cloud mass (md^o)- The regions with nicofl and 
mci.o are then regarded as the core and cloud, respectively. The threshold densities (pth.corc and 
Pth,cl) vary with time, and are obtained from eq. (9) by means of a bisection method. Unless noted 
otherwise, we will choose rrid as the threshold mass to obtain the global quantities. 

To measure the evolution of the cloud shape, we define the eff^ective radii normal to and along 
the cylinder axis, a and c, respectively, as 



5 \ V2 



1/2 



where 



m. 



m. 



'p>Pth 
-1 / „„2 



pzdV, {z^) 



m. 



pz^dV, 



p>Pth 



th 



pw^dV . 



(10) 

(11) 
(12) 



p>Pth 



To measure the cloud destruction and turbulence, we define the velocity dispersions in the 
axial and radial directions, respectively, as 



2 \l/2 



where 



Svz = (.{vl) - {v^)Y' and Sv^ = {vl) 



{vl) = m^^ / pv^dV , (vl) = m^^ / pv^dV 
•Jp>Pth •Jp>Pth 



The mean velocity, (vz), is defined as 



(vz) = / pVzdV 



(13) 
(14) 

(15) 



where Vz is measured in the frame of the unshocked cloud. Finally, the mean density, (p), is defined 



as 



(p) = Kh' / pdy ' 



(16) 



p>Pth 



where Vth is the volume of the region with p > pt^. In the following, the initial mean density is 
refered to as (po) and the mean initial density contrast is defined as (x) = (po) / Pio- 
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3. 



Numerical Simulations 



We solve the two-dimensional Euler equations of hydrodynamics using an Eulerian adaptive 
mesh refinement (AMR) hydrodynamic code with a second-order Godunov, dimensionally-split 
approach (Truelove et al. 1998; Klein 1999). In this code, as a computation proceeds and one 

wishes to resolve certain features within the flow, fine grids arc dynamically created at locations 
meeting a pre-determincd criterion, initialized from interpolated existing coarse grid data, and 
evolved. This entire process is continually repeated, allowing us to create and destroy fine grids in 
a general and controllable fashion, wherever and whenever the need arises (Truelove et al. 1998). 

In this paper, the entire computational domain is covered with a base grid (level 1) and we use 
2 additional levels of refinement (levels 2 and 3). The entire computation domain is covered with 
a base grid (level 1). Finer grids (levels 2 and 3) are created or destroyed by using two refinement 
criteria, a density criterion and density gradient criterion. The density criterion leads to refinement 
of cells in the base grid to level 2 when the density in these cells exceeds a threshold density, p^ef- 
In all the calculations presented in this paper, wc set prcf = 1-2 pn, where pn is the density of the 
post-shock ambient gas. This guarantees that the initial cloud (po > 2 pio) is refined to at least 
level 2 after passage of the shock. We followed several models with different p^ef and found that 
the numerical results are essentially identical for p^ef ^ l-5pii (i.e., the global quantities defined 
in the previous subsection are converged to within a few percent). The density gradient criterion 
leads to refinement of levels 1 or 2 when the local density scale height normalized by the cell size 
becomes smaller than a threshold value, Zth; 



where Axj denotes the grid spacing of the i-th grid and the density gradient dp/dx\j is computed 
at the cell center. We typically adopt l^^ = 20 in this paper because the global quantities are 
converged to within a few percent for Ith > 10. This criterion guarantees that shock fronts and 
contact discontinuities are refined to level 2 in the region of p < pref and to level 3 in the region of 

P > Pref ■ 

Our grids are Cartesian and properly nested; i.e., a fine grid must be surrounded by at least 
one layer of coarser cells. The only exception occurs at domain boundaries, where fine and coarse 
interfaces may coincide. The increase in resolution between levels, referred to as the refinement 
ratio, is taken to be 4 for every level. 

The spatial resolution of the calculation is specified by the number of finest grid points that 
covers an initial core radius Vco- We use 120 finest grid points per core radius Tco as our standard 
resolution because we have demonstrated that with this resolution the evolution of the global 
quantities is well-converged. We refer to this resolution as Ri2o- More generally, we define the 
resolution of a calculation hy Rn, where the subscript N is the number of cells per Tco at the finest 
level of refinement. The spatial resolutions we used in this paper and several characteristics of the 
initial clouds are summarized in Table 1. 




(dp/dx\j)Axj 



(17) 
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Since our model cloud does not have a sharp boundary between the cloud and intercloud 
medium, we do not explicitly separate the cloud and intercloud medium as different fluids; i.e., we 
use a single-fluid version of the AMR code in which the ratio of the specific heats, 7, is always 
identical for both the cloud and intercloud medium. By contrast, KMC used a two-fluid code that 
treated the cloud and the intercloud medium as different fluids. When we approximate the effect of 
radiative losses by setting 7 ~ 1, we must assume that both the cloud and the intercloud medium 
are radiative, whereas KMC were able to consider a radiative cloud in an adiabatic intercloud 
medium. As a result, our results for the radiative case differ from theirs (see §4.5). 



3.1. Convergence Tests 

In hydrodynamic simulations, it is very important to demonstrate that the calculations are 
performed at spatial resolutions that resolve important features of the hydrodynamic evolution such 
as the growth of hydrodynamic instabilities (see also KMC and Klein & Woods 1998). Comparing 
the time evolution of the global quantities at different spatial resolutions, KMC showed that at 
least 10^ cells per cloud radius are needed to follow the interaction of shock waves with uniform 
clouds accurately. Here, we carry out a similar study for smooth clouds. 

To perform the convergence study, we calculate two models for the shock-cloud interaction 
with {n,x,M) = (2,10,10) and (8,10,10) with different spatial resolutions (-R30) -^60; R120, R240, 
i?480, and i?960 for n = 8 and Reo, R12O; cind i?48o for n = 2) and compare the global 

quantities. Figure 1 shows how the errors at the end of the calculation (t = 3 tec for n = 8 and 
t = 6tcc for n = 2) depend on the spatial resolution. The errors at a resolution Rj\j arc defined as 
the differences between the global quantities at Rn and those at the finest resolution, Rf. In other 
words, for a global quantity Qjv measured at Rn, the error is 

Ari _ \Qn - Qf\ .-.^s 
AQ^ - i^^i • (18) 

The convergence test indicates that the errors of the velocity dispersions generally decrease in 
proportion to the inverse of the cell number N, which implies first-order accuracy. On the other 
hand, the convergence of the other quantities is better than that of the velocity dispersions; e.g., the 
errors in the radii a and c scale as iV~^. The slow convergence of the velocity dispersions may be 
related to the fact that the code accuracy is reduced to first order at discontinuities such as shocks 
and contact discontinuities, where substantial velocity differences are generated. All the quantities 
converge to within about 10 % for the spatial resolutions of -R120 or finer. Hence, we conclude that 
at least ~ 10^ cells per rco are necessary for a converged calculation. In the following, we adopt 
i?i2o as the standard spatial resulution. 
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4. Numerical Results 

In this section, we describe the overaU evolution of clouds with smooth density profiles. In 
§4.1 and §4.2, we consider both clouds with a relatively steep density gradient (n = 8) and with 
a shallow density gradient (n = 2) as representative cases of nonradiative (7 = 5/3) clouds with 
smooth density profiles. In these models, the other parameters are fixed at a density contrast of 
X = 10 and a Mach number M = 10. The dependence of the results on the parameters n and x is 
summarized in §4.3. In §4.4 we discuss the applicability of Mach scaling (see §4 of KMC) to smooth 
clouds. The effect of 7 is studied in §4.5. Finally, the evolution of cylindrical clouds is compared 
with that of spherical clouds in §4.6. 

In describing our results, we shall unit of length, ^ & unit of density, and tec 

as a unit of time. Thus, our model is specified by four dimensionless parameters: n, X: M, and 7. 
The evolution time is measured from the epoch at which the initial shock has reached to the point 
z = — rco on the z-axis. In Table 1, we summarize the parameters of the models presented in this 
paper. 

According to previous studies, the overall evolution of shocked uniform clouds can be divided 
into four different stages: (1) the initial transient stage, in which an incident shock is transmitted 
into the cloud, and at the same time a bow shock or bow wave propagates into an upstream region; 
(2) the shock compression stage, in which the cloud is compressed mainly in the z-direction; (3) 
the postshock expansion stage, in which the shock-compressed cloud expands laterally, and (4) the 
cloud destruction stage, in which the shocked cloud is destroyed primarily by Kelvin-Helmholtz 
and Rayleigh- Taylor instabilities. As shown below, while the overall evolution of smooth clouds 
can also be divided into these four stages, there are significant differences in the morphology and 
evolution. 



4.1. A Cloud With A Steep Density Gradient: ASS {n = 8,x = 10, and M = 10) 

To begin with, we consider a cloud with n = 8, x = 10) ^-iid 7 = 5/3 interacting with a 
Mach 10 shock. In this model, about 70 % of the cloud mass is concentrated in the central core 
{fnco,o ~ 0.7mci,o)- There are several significant differences between the evolution of a uniform 
shocked cloud and a smooth shocked cloud, as we now describe. 

In Figure 2, we present snapshots of the density distributions at seven different times. At the 
first encounter with the incident shock, the cloud is compressed mainly in the z-direction (the initial 
transient stage). In this stage, a shock reflects from the front of the cloud, propagating away from 
the cloud into the upstream postshock flow. The reflected shock develops into a bow shock, which 
can be recognized as a density jump in the upstream flow in Figures 2b and 2c. The bow shock 
eventually becomes a weak acoustic wave at later stages of evolution (Spitzer 1982). In fact, the 
bow shock is not so prominent in Fig. 2d as in Figures 2b and 2c. In Figure 2b, the shock passing 
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over the cloud envelope converges on the z-axis behind the cloud (at zjrco ^1.5 — 2 in Fig. 2b), 
reflects off the z-axis, and by a double Mach reflection drives a shock back into the shocked cloud. 
Thereafter, the shock transmitted into the front of the cloud (the cloud shock) interacts with the 
reflected shock after the cloud shock has passed through the high density region. (For the uniform 
case, this process occurs inside the cloud. Compare Fig. 2b and Fig. 2 of Klein et al. (1990).) This 
reflected shock compresses the shocked cloud from the rear in the z-direction (the shock compression 
stage). At this stage, the shocked cloud is deformed into an arc- like structure without any sign 
of the growth of instabilities (Fig. 2c). The cloud then expands laterally because the pressure is 
lower at the sides of the cloud compared to that on the axis (the postshock expansion stage, see also 
Nittman et al. (1982)). At the time shown in Figure 2d, the Kelvin-Helmholtz instability begins 
to develop along the slip surface, and eventually destroys the shocked cloud (the cloud destruction 
stage, see Figs. 2e, 2f, and 2g). The cloud has broken up into small fragments by the time shown 
in Figure 2f. 

Since the velocity of the cloud shock is smaller in the higher density region (since v oc p"^/^), a 
significant velocity gradient develops in the cloud envelope. In the shock compression and postshock 
expansion stages, this velocity gradient in the cloud envelope steepens the initial smooth density 
profile at the front of the cloud, resulting in a sharp density jump (Fig. 2b). We refer to this density 
jump as a "s/ip surface" because the velocity is sheared parallel to its surface. The evolution of the 
slip surface is crucial in understanding the evolution of smooth clouds, because the growth rate of 
the Kelvin-Helmholtz instability across the surface depends on the thickness of the surface itself. 

The convergence of the intercloud shock on the ^;-axis creates another strong reflected shock, 
which propagates downstream. It interacts with the intercloud shock to produce a Mach-reflected 
shock, which develops into a double Mach-reflected shock with two triple points. A powerful 

supersonic vortex ring forms just behind the Mach-reflected shock, and is carried away from 
the cloud with almost the same speed as the incident shock. The vortex ring propagates up to 
{w/rco, z/rco) — (0.2,5) in Figure 2c, and moves up to {w/rcojz/rco) — (0.5,11.4) and (0.6,17.6) 
in Figures. 2d and 2e, respectively. 

Immediately after the shock passage, the morphology of the shocked cloud is significantly 
different from the uniform case, for which the lateral expansion occurs more rapidly and the Kelvin- 
Helmholtz instability starts to grow immediately after the shock passage. For example, for the 
uniform cloud, the ripple generated by the instability already apprears along the cloud surface 

at t ~ Itcc [see Figs. 1- 4 of Klein et al. (1990)]. and the arm-like structure starts to break up 
away from the cloud by t ~ 3 tec- For the smooth cloud, the cloud is deformed into an arch-like 
structure due to the shock compression, without any signs of the Kelvin-Helmholtz instability even 
at 2 tec- The retarded growth of the instability is due to the fact that the instability grows faster for 
shorter wavelengths (the growth rate oc k~^, where k is the wavenumber), but is strongly suppressed 
for wavelengths shorter than the thickness of the shear layer [See Chapter 11 of Chandrasekhar 
(1961)]. Therefore, cloud destruction is retarded until a density discontinuity develops; in this 
case, the density discontinuity is associated with the slip surface. Also, the Richtmyer-Meshkov 
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instability, which grows hnearly in time for small amplitudes, as opposed to the exponential growth 
of the Kelvin-Helmholtz and Rayleigh- Taylor instabilities in the linear regime, does not develop 
appreciably on the z-axis (Compare Fig. 2b and Fig. 5 of Klein et al. (1990)). In §6 and 7 we shall 
discuss in more detail the formation of the slip surface and the subsequent destruction of the cloud 
by the Kelvin-Helmholtz instability. 

4.2. A Cloud With A Shallow Density Gradient: AS2 (n = 2, x = 10, and M = 10) 

Next, we consider a very smooth cloud with n = 2. The other parameters are the same as 
those of the previous model: x = 10) 7 = 5/3, and M = 10. The density profile of this cloud 
resembles that of a Bonner-Ebert sphere, or isothermal equilibrium sphere, which has a density 
profile p oc in its envelope. The cloud mass is less centrally concentrated than in the previous 
case (mco,o = 0.10mci^o)j only about 10 % of the total cloud mass is in the central core. Therefore, 
the smooth density profile has a greater effect for this model than for the n = 8 model presented 
in §4.1. 

The shallow density gradient makes the interaction of the shock with the cloud milder than 
in the previous case. In Figure 3 we show snapshots of the density distributions of the shocked 
cloud. At the initial transient stage, a bow wave forms in the upstream flow, as can be seen at 
z/rco — — 2 ^ — 1 in Figs. 3b and 3c. However, it does not evolve into a shock, in contrast to the 
previous case and to the uniform cloud case. The reflected shock in the downstream flow is also 
weak. Because the effective density contrast between the cloud and intercloud medium is small 
((Pci) / Pio ~ 2.9) due to the small mass concentration, the mean relative velocity between the cloud 
and intercloud medium is also small. As a result, the intercloud shock does not converge on the 
z-axis at the rear of the cloud but instead interacts directly with the cloud shock (compare Figs. 
2b and 3b). Thus, shock compression from the downstream side is weak, resulting in a slow lateral 
expansion. The interaction also produces a single Mach-reflccted shock with a single triple point, 
which propagates away into the downstream flow from the cloud with almost the same speed as 
the intercloud shock (Figs. 3c, 3d, and 3e). Although the intercloud shock does not reflect from 
the z-axis, the postshock flow tends to converge on the z-axis, producing a vortex ring near the 
z-axis. This vortex ring propagates with a speed slower than the Mach reflected shock by a factor 
of ~ 0.75. 

Another effect of the shallow density gradient is a significant retardation in the formation of 
the slip surface. Because of the smaller mean relative velocity between the cloud and intercloud 
medium, it takes more time to form the slip surface. The Kelvin-Helmholtz instability grows very 
slowly due to both the small velocity shear and the retardation of the slip surface formation. An 
undulation along the slip surface does not appear until t ~ 9icc (Fig- 3e), whereas in the n = 8 
case, it alreadly appears at i ~ 4 tec (Fig- 2d). Thus, it takes much more time for the cloud to be 
completely destroyed for n = 2 than for n = 8. 
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4.3. Dependence on n and % 



As mentioned in §2, our model is specified by four parameters: tlie power index of the density 
profile (n), the density contrast (x), the strength of the incident shock (M), and the ratio of specific 
heats (7). Here, we consider the effects of varying n and %, which determine the mass distribution 
of the initial cloud. The dependence on M and 7 will be discussed in subsequent subsections. 

First, we consider the effect of the density gradient at the cloud boundary (see Figs. 2 and 
3 and Tables 2 and 3). For smaller n, the interaction of the shock with the cloud is milder. For 
example, during the initial transient stage, the bow shock is weaker for smaller n, and the reflected 
shock and supersonic vortex ring formed at the downstream are less powerful. This is in part due 
to the fact that the mean relative velocity between the cloud and the intercloud medium decreases 
with n. For smaller n, the initial mean density contrast between the cloud and the intercloud 
medium ((x) = (pd) /pio) is smaller (see the seventh column of Table 1), and thus the relative 
velocity just after the initial shock passage is smaller, since 



for strong shocks. This in turn delays the formation of the slip surface and therefore the onset of 
the Kelvin-Hclmholtz and Ray Icigh- Taylor instabilities. As a result, it takes longer for clouds with 
smaller n to be destroyed: The cloud destruction times ijcst for the n = 2 models are longer than 
for the n = 8 and n = 00 models by factors of ~ 2 and ^ 3, respectively (see §7 for the definition 
of idest)- 

The radial expansion is milder for smoother clouds; the radial expansion speed decreases with n 
because the relative velocity between the cloud and intercloud medium decreases with n (see §5.2). 
Also, the maximum effective radius normalized by the initial radius, a/oo, decreases monotonically 
with n, reflecting the fact that the shock-cloud interaction is milder for smaller n. 

Next, we consider the effect of varying the density contrast x- Figure 4 compares the snapshots 
of the X = 100 clouds with n = 8 and n = 2, and Table 3 summarizes the global quantities for 
X = 100 over a range of n. The most important effect of increasing the density contrast x is that 
it leads to an increase in the relative velocity between the cloud and the intercloud medium (vj-ei 
increases with (x) according to eq. 19 and (x) increases with x according to Table 1). This leads 
to a longer drag time (see eq. [2.6] of KMC). More significantly, because the Kelvin-Helmholtz and 
Rayleigh- Taylor instabilities both depend on the relative velocity between the cloud and intercloud 
medium, (for the Rayleigh- Taylor instability, this is because the higher relative velocity leads to a 
greater acceleration) the higher relative velocities at higher x lead to faster growth of the instabilities 
there. For example, for n = 2, a small-scale undulation of the slip surface due to the Kelvin- 
Helmholtz instability already appears at t ~ 4 fee for x = 100, while it does not appear until 
i ~ 8tcc for X = 10 (compare Figs. 4b and 3d). At x = 100, the Rayleigh-Taylor instability leads 
to substantial destruction of the cloud near the z-axis, whereas the Rayleigh-Taylor instability 
does have significant effects at x = 10 (see Figs. 2 and 3). Our simulations also show that the 
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cloud destruction timescale measured in units of cloud-crushing times depends only weakly on 
For the x = 10 and 100 clouds with the same values of n, the results for tdcst /^cc coincide within 
errors of ~ 10 percent (see Tables 2 and 3). This is consistent with the analytic estimate of KMC, 
which suggests that the growth times of the Kelvin-Helmholtz and Rayleigh- Taylor instabilities are 
comparable to the cloud-crushing time (see §2 of KMC). 

Another important effect of increasing the density contrast is that the intercloud shock is able 
to converge onto the ^;-axis for very smooth clouds (n = 2) at x = 100, whereas it cannot do so at 
X = 10. As a result, at x = 100 a strong reverse shock is driven back into the cloud from the rear, 
compressing the cloud in the axial direction. After the reflected shock collides with the transmitted 
shock in the cloud, the cloud undergoes a steady expansion in the axial direction that increases 
with X- At later times, the Rayleigh- Taylor instability also contributes to the larger effective radius 
in the axial direction. On the other hand, the radial expansion seems to reach an asymptotic value 
(see §5). 



4.4. Dependence on M 

4- 4-1- Strong Shocks: Mach Scaling 

For strong shocks (M — > oo), the hydrodynamic equations for an inviscid gas are invariant 
under the transformation 

t-^tM, v/M, P P/M"^, (20) 

with the position and density unchanged. This is referred to as "Mach scaling". For uniform 
clouds, KMC demonstrated that Mach scaling is reasonably valid. Here, we investigate whether 
Mach scaling is applicable to smooth clouds. 

In Figure 5, we compare the density distributions of a cloud with n = 8, x = 10 and 7 = 5/3 
for two different Mach numbers, M = 100 and M = 1000. The clouds are depicted at two different 
times: t = 3.5tcc (Figs- 5a and 5b) and t = GAtcc (Figs. 5c and 5d). At the first epoch (Figs. 5a 
and 5b), the cloud morphologies are essentially the same. At the second epoch (Figs. 5c and 5d), 
morphological differences begin to appear at small scales, even though the overall morphologies 
remain similar. Those morphological differences become more significant at later stages. These 
deviations most likely originate from truncation errors and the artificial viscosity, which do not 
scale as equation (20). 

In spite of the deviations in the cloud morphology, the global quantities agree very well even 
at late times if all the variables are scaled as in equation (20). For example, for the M = 100 and 
M = 1000 cases, the effective cloud radii and velocity dispersions coincide within about 5 percent 
and 10 percent, respectively, until the final stage of the computation {t ~ 10 tec)- (Values of these 
quantities at t = t ui ~ 5 tec Q-rG listGci in Table 4.) TI16 clia-ractcristic tiinescalGs such, os ^clrag? ^dest? 
^mix, and tm are in excellent agreement within errors of a few % for strong shocks. 
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4- 4- 2- Weak Shocks 

Here we consider the weak shock case, in which Mach scaUng is not applicable. We follow two 
cases with initial parameters that are the same as those of §4.1 and §4.2, except that the Mach 
number is M = 1.5. 

In Table 4, we summarize the global quantities of the models with (n, x, M) = (8, 10, 1.5) and 
(2, 10, 1.5). Our numerical simulations show that the overall evolution of the weak shock case is 
qualitatively similar to that of the strong shock case. An important difference between the weak 

and strong shock cases is that cloud drag is less efficient in the weak-shock case, in part due to 
slower lateral expansion of the shocked cloud. For example, for n = 8 and 2, the drag times are 
17.2 tec and 21.4 tec, respectively, for the weak-shock case, whereas they are 3.3 tec «md 7.5 ico 
respectively, for the strong-shock case. 

Furthermore, the relative velocity between the cloud and intercloud medium is smaller in the 
weak-shock case, even after allowing for Mach scaling. This smaller relative velocity leads to slower 
growth of the Kelvin-Helmholtz instability (and thus slower cloud destruction). For example, for 
n = 8, the ripple along the slip surface due to the Kelvin-Helmholtz instability begins to appear at 
t ~ 12 tec in the M = 1.5 case, while it already appears at t ~ Atcc for M = 10. Thus, the cloud 
destruction time is about 2 — 3 times longer in the weak-shock case. 

The supersonic vortex ring is weak or absent for weak shocks. For M = 1.5 wc did not observe 
the formation of a supersonic vortex ring for either n = 8 or n = 2, although the interaction between 
the reflected shock and the incident shock did generate a single triple point. Also, for the weak 
shock, a bow wave forms instead of a bow shock. Despite these qualitative differences between the 
strong and weak shock cases, the overall evolution of the weak shock case is quite similar to that 
of the strong shock case. We thus conclude that varying the strength of the incident shock does 
not alter our main conclusions on the shock-cloud interaction. 

4.5. Dependence on 7: Effects of Radiative Losses 

In the above subsections, we considered the evolution of nonradiative clouds, in which 7 was 
set to 5/3. Here, we consider the effects of radiative losses, which we model by setting 7 = 1.1. 
Such a model is applicable for relatively large HI clouds and molecular clouds (e.g., KMC, Larson 
2005) We recalculate the evolution of two models with (n,x,M) = (8, 10, 10) and (2, 10, 10). The 
numerical results are summarized in Table 5. Note that in our calculations, 7 is 1.1 for both the 
cloud and the intercloud gas. In contrast, KMC set 7 = 5/3 for the intercloud gas in their radiative 
calculations. 

The softer equation of state leads to substantially greater compression; e.g., for the radiative 
model with n = 8, the effective radii, a and c, are reduced by as much as factors of 2.3 and 2.6, 
respectively, whereas for 7 = 5/3, a and c are reduced by factors of 1.3 and 1.8, respectively. Thus, 
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the cross section of the radiative cloud is smaller than that of the nonradiative cloud, which lessens 
the efficiency of the cloud drag by the ambient flow. On the other hand, the post-shock intercloud 
density is higher than in the nonradiative case by a factor of about 4.7 (for M = 10). This higher 
intercloud density gives a greater drag force on the shocked cloud than in the nonradiative case (see 
eq. 21 below). In our calculations, the latter effect is more important, and as a result, the radiative 
drag times are much shorter than the nonradiative ones for both the n = 8 and n = 2 models. 
Because of this greater acceleration by the ambient gas, the relative velocity at the cloud boundary 
is rapidly reduced, leading to milder cloud destruction by the Kelvin-Helmholtz and Rayleigh- Taylor 
instabilities (see Table 5). The higher compression in the radiative case may facilitate gravitational 
collapse in the shocked cloud (§9). 



4.6. Effect of Cloud Geometry 

To investigate the effect of the cloud geometry, we consider cylindrical clouds. We adopt a 
Cartesian coordinate system {x, y, z) and assume that the physical quantities are uniform along the 
z-axis, which is identical to the cylinder axis. Mirror symmetry is assumed with respect to the x-z 
plane. The shock propagates in the y-direction. Hence, the x- and z-axes in this case correspond to 
the zu axis in the spherical case, and the y-axis in this case corresponds to the z-axis in the spherical 
case. In Table 6, we summarize the global quantities of the cylindrical clouds with (n, x, M) = (2, 
10, 10) and (8, 10, 10). Comparison between the cylindrical and spherical clouds indicates that the 
cloud shape becomes more flattened and longer in the transverse direction for the cylindrical cloud. 
For example, for the n = 8 case, the effective radius in the transverse direction reaches a/rco — 2.5 
in the cylindrical cloud, which is about 20 % larger than that of the spherical cloud with the same 
n. The somewhat greater expansion of the cylindrical cloud causes more efficient drag. Therefore, 
the cloud is accelerated more rapidly for the cylindrical cloud, leading to a slightly slower cloud 
destruction due to smaller relative velocities between the cloud and intercloud medium. This slower 
cloud destruction also causes a slower increase in the axial cloud size for the cylindrical cloud; at 
t = 12 tec, it reaches about c/vco ~ 4.0 in the cylindrical cloud, while it reaches about c/vco ~ 4.7 in 
the spherical cloud. In spite of these quantitative differences, the overall evolution of the cylindrical 
clouds is quite similar to that of the spherical clouds. Thus, we conclude that varying the cloud 
geometry does not alter our main conclusions obtained from the spherical calculations. 



5. Cloud Drag 

As shown above, the slip surface is subject to the Kelvin-Helmholtz and Rayleigh- Taylor 
instabilities, and as a result, the entire cloud is eventually destroyed. The acceleration of the cloud 
is responsible for the growth of the Rayleigh-Taylor instability, and the time history of the velocity, 
which is determined by the acceleration, governs the evolution of the Kelvin-Helmholz instability. 
Therefore, it is very important to quantitatively understand the cloud acceleration. For uniform 
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clouds, KMC constructed a simple analytic formula for the cloud velocity that reproduced the 
numerical results quite well. Here, we extend their analytic model to smooth clouds. 

5.1. Analytic Formula for the Cloud Velocity 

According to §2 of KMC, the equation of motion for a shocked cloud can be expressed as 

mcl,0^ = -\cDP^lvlA{t) , (21) 

where mdfi is the cloud mass, Vc is the mean velocity of the cloud relative to the velocity of the 
shocked intcrcloud medium, pn is the density of the post-shock ambient medium, C/j is the cloud 
drag coefficient, and A{t?) is the time-dependent cross section of the cloud (sec Figs. 2 and 3). As 
shown below in Figure 6, the cloud radius (a) is reduced slightly in the shock compression stage, 
and then increases significantly with time during the post-shock expansion stage. Thereafter, it 
does not change much after it reaches its maximum. Here, we approximate the evolution of the 
cross section of the cloud as 

A{t) = TTO^ = TT [al + {Cctf] t < tm (22) 
= Tr[al + [Cctrr^f] t>tm, (23) 

where ao is the effective value of the initial cloud radius and Cc is the expansion velocity in the 
tZ7-direction. In this equation, we assume that the cloud expands at a constant speed, Cc, until it 
reaches its maximum transverse extent. We define the maximum expansion time, tm, to be the time 
at which the effective radius in the oj-direction (a) has increased to 90 % of its maximum value. 
The maximum expansion times obtained from the numerical calculations are listed in Tables 2 — 
6. Equation (23) is identical in form to the expression given by KMC, but there is an important 
difference: Cc is not necessarily equal to the sound speed in the cloud, but instead depends on 
model parameters. 

In order to obtain a result that is directly comparable to that of KMC, we work in the frame 
of the shocked intcrcloud medium, denoted by a prime. Integrating equation (21) with respect to 
time, we obtain the following form of the cloud velocity, 

(t < tm) (24) 
{t > tm) (25) 
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We have assumed that the shock is strong (M 1) and made the approximation Tco — Vgtcc, 
where Vs — vi^/x^^"^ is the velocity of the transmitted cloud shock. This result is identical to that 
of KMC except that the drag coefficient C^i is replaced by an effective drag coefficient that 
can be much less than Co for very smooth clouds (see below eq. 4). In the case considered by 
KMC (a uniform cloud [TZ = 1] in an intercloud medium with 7j = |), = Cd- The initial cloud 
velocity relative to the shocked ambient gas is v'^q = —2vb/{'y + 1). In the following, we set the drag 
coefficient Cd = 1 (see KMC). 



5.2. Comparison with Numerical Simulations 
5.2.1. Cloud Shape 

In the left panels of Fig. 6, we show the time evolution of the rms cloud radii, a and c. The 
cloud radius in the transverse direction, a, decreases with time in the shock compression stage, and 
then increases in the post-shock expansion stage, finally reaching a maximum. Thereafter, it shrinks 
by about 5%-30%, depending on n and X) with the greater shrinkage occurring for the larger values 
of X- On the other hand, the axial radius, c, keeps increasing with time. Thus, the shocked cloud 
becomes elongated in the direction of shock propagation. The initial reduction in the transverse 
radius due to shock compression happens at a later time for smaller n; for (x, M) = (10,10), the 
radius a takes its minimum value at t ~ 1.7 tec for n = 8 and 3.2 tec for n = 2. As n and x increase, 
the cloud expands more rapidly, the aspect ratio of the shocked cloud, c/a, increases more rapidly, 
and the maximum value of a/ao increases. 

For a uniform cloud, the expansion speed of the cloud can be approximated as the sound speed 
of the shocked cloud material (KMC), 



Cr. Cc 



X(7+1)^J ■ ^ ^ 



For X = 10) this gives Cc = 0.56^5 for 7 = 5/3. For smooth clouds, the expansion speeds are 
smaller: C^Vs ~ 0.07 , 0.23, and 0.25 for n = 2, 4, and 8 with x = 10) respectively. The increase 
of expansion speed with n is due to the increase in relative velocity between the cloud and the 
intercloud medium as the cloud boundary becomes sharper and the mean density contrast larger; 
this in turn leads to a larger pressure difference between the shocked cloud and the intercloud 
medium. 

The expansion speed also tends to increase with x- 1^ addition to the effect just described, 
the shock driven into the rear of the cloud augments the pressure in the cloud interior due to the 
transmitted shock. As a result, the radial expansion due to the pressure difference between the 
cloud and the side of the cloud is somewhat faster. Indeed, for n > 8 and % S> 10, the expansion 
speed Cc is approximately equal to the adiabatic sound speed in the cloud, just as in the case of 
uniform clouds. 
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5.2.2. Cloud Velocity 

In the right panels of Fig. 6, we compare the analytic result for the cloud velocity given in 
equation (25) with the value obtained from the numerical calculations. As shown in the figures, 
the expansion velocity Cc depends on n and x- Since we determine Cc by fitting the evolution 
of the cIoTid radius a in the time interval t < tm, our result for the cloud velocity is really only 
semi-analytic. 

For clouds with a large density gradient, i.e., n = 8, our analytic model of the cloud velocity 
agrees well with the numerical results for both x = 10 and 100. For x = 10; th^ analytic cloud 
velocities are within about 25 % and 10 % of the numerical values for t < tdrag a-^d t > tdrag) 
respectively. The cloud velocity of the n = 8 cloud drops more slowly than the uniform case 
because its radial expansion is slower, and thus its cross section is smaller. For x = 100, the 
numerical values and analytic model agree to within 10 % until the radius a reaches its maximum. 
Thereafter, the deviation increases as the cloud undergoes a brief bout of acceleration, presumably 
due to the effects of the reverse shock. A similar effect is observed in some other models [see e.g., 
panel (d)], where it leads to a deviation from our model that is only temporary. 

For smooth clouds with ra = 2, the initial reduction in the effective cloud radius associated with 
the cloud compression significantly affects the cloud acceleration. Therefore, for such clouds we 
replace oq in eq. (23) with the time-averaged radius over the time interval to < t < ti, where to is 
the time at which the shock first reaches the cloud surface on the 2;-axis, and ti is the time at which 
the cloud radius expands back to its initial value. This procedure reduces ao in eq. (23) by about 
10% and 17% for x = 10 and 100, respectively. Even with this adjustment, the agreement between 
the analytic model and the numerical simulations at early times is not as good as for the n = 8 
cases. This is because a significant fraction of the mass of the cloud is not much denser than the 
intercloud medium (the mean density contrast between the cloud and the intercloud medium is only 
(p) / Pio ^ ^ for the n = 2 case), and as a result the transmitted shock gives this gas a significant 
velocity. However, in deriving equation (25), we assumed that the cloud velocity is essentially zero 
immediately after the initial shock passage. At later times, the agreement between the analytic 
model and the numerical solution becomes much better in the n = 2 case: for 5tcc ^t < 20 tec j the 
agreement is within about 20 % for both x = 10 and 100. Note that the cloud was still expanding 
at the end of the calculation for x = 100 (t < tm), so we were unable to fully test the model in this 
case. 

In summary, our analytic model of the cloud drag can reproduce the numerical results reason- 
ably well over a wide range of initial conditions. 
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6. Formation of the Slip Surface 

As shown in §4, smooth clouds are completely destroyed by Kelvin-Helmholtz and Rayleigh- 
Taylor instabilities, whose growth rates depend on the wavelength of the perturbation. According 
to linear theory, the growth timescales of these instabilities are estimated as t-Kn/tcc ~ vi,/{v,:eikrco) 
and tKi/f-cc = l/(^'"co)~^^^; respectively (see §2 of KMC), where k is the wavcnumber of the 
perturbation. In both cases, the instability grows faster on smaller scales. Furthermore, linear 
theory shows that perturbations with wavelengths smaller than the thickness of the shear layer are 
stabilized significantly. Thus, the growth of the instabilities is closely related to the development 
of the slip surface. In this section, we investigate how a slip surface develops in a smooth cloud. 

In Figures 7a — 7c, we represent the distributions of p, Vr, and along r = 0.5rco at five 
different times for the cloud with small density gradient, i.e., n = 8 and x = 10- When the trans- 
mitted shock passes through the cloud (dash-dot-dot-dot curves), the maximum density increases 
by a factor of ~ 4, as inferred from the jump condition, and then decreases as the cloud expands 
laterally. At the first encounter of the shock, the velocity shear layer is formed at the cloud enve- 
lope with a thickness /ibd- As can be seen in panel (c), the cloud envelope at the upstream side 
moves faster than the high density core, catching up with the high density core. As a result, the 
radial velocity profile bcomes steep there [panel (b)], and thus the slip surface forms. At the same 
time, the density profile becomes steep and therefore the slip surface is also identical to a contact 
discontinuity. Here, we define the formation time of the slip surface, tglip) as the time at which the 
radial velocity gradient reaches a first local maximum at the upstream side of the shocked cloud. 
In this model, tglip is estimated as tgUp ~ 1.71 tec- At r = 0.5rco, where we measured the physical 
quantities, the slip surface is almost perpendicular to the z-axis, and thus the radial velocity dif- 
ference across the slip surface can be regarded as the total velocity difference there. We note that 
the definition of tsUp depends on the position where the radial velocity is measured. However, as 
long as we measure the radial velocity gradient in the range of 0.3 < r/vco ^ 0.7, we get almost 
the Sciine value. Thus, we mecisurc tslip t — O.^t^q for all models. At t ~ 3tcc, the reverse shock 
compresses the cloud from the rear, and the cloud is squeezed from both sides by the reverse shock 
and the post-shock intercloud flow. Thus, the density gradient of the slip surface becomes very 
steep as can be seen in the solid curve in panel (a). The slip surface forms close to the core, and it 
can form only after the transmitted shock has reached the core. 

For the cloud with small density gradient, i.e., (n, x) = (2, 10), the shear layer has a wider 
thickness than for n = 8, at the first encounter of the shock. The lateral expansion of the core is a 

little faster for the smoother cloud, which leads to a smaller velocity difference between the core and 
intercloud medium. Because of this smaller velocity difference, the formation of the slip surface is 
retarded for the smooth cloud. In this model, tgHp is evaluated as t^np = 4.92 tec- Since the reverse 
shock is weak in this model, the additional compression due to the reverse shock observed in the 
n = 8 case does not happen. 

When the slip surface forms, the Kelvin-Helmholtz instability begins to develop and produce 
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ripples along it. For all the models calculated in this paper, the velocity dispersions in the axial 

and radial directions, Sv^ and 5vz, begin to rise at the epoch of ripple formation due to the Kelvin- 
Helmholtz instability. The time evolution of the velocity dispersions is shown in Figure 8 for two 
models with (n, x, M) = (2, 10, 10) and (8, 10, 10). For the n = 8 case, the velocity dispersions in the 
w- and z-directions begin to increase at t c± 2.3 tec 3.2 tec, respectively. For the n = 2 case, they 
begin to increase at t ~ 7 tec and 8 tec, respectively. These times are in good agreement with the 
times of ripple formation (3.5 tec and 8.4 tec for n = 8 and 2, respectively). This indicates that the 
small-scale perturbations due to the Kclvin-Helmholtz instability indeed generate random motions 
of the small fragments. We note that the velocity dispersions in the axial and radial directions reach 
about 10 % of the incident shock velocity over a wide range of initial model parameters. They tend 
to be slightly larger for clouds with steeper density profiles (larger n — see Tables 2 through 6). 

7. Cloud Destruction 

The transmitted shock deforms a smooth cloud into an arc-like structure bounded by a slip 
surface. The slip surface is subject to the Kelvin-Helmholtz and Rayleigh-Taylor instabilites, which 

produce fragments that arc swept away by the flow. Although their growth rates arc smaller, 
fluctuations with larger wavelengths break the cloud into relatively large fragments, destroying the 
cloud more efficiently. On the other hand, fluctuations with smaller wavelengths contribute more 
efficiently to mixing the cloud material with the ambient gas. To evaluate these two processes of 
cloud destruction, we consider two different timescales: the cloud destruction time (tdest) f'^d the 
mixing time (tmix), which are defined in §7.1 and §7.2, respectively. 

7.1. Cloud Destruction Time: tdest 

When fluctuations with wavelengths comparable to the cloud size grow in time, they break the 
cloud into relatively large fragments. This large-scale fragmentation decreases the mass of the main 
body of the shocked cloud efficiently. To make a quantitative estimate of large-scale fragmentation, 
we follow KMC and define the cloud destruction time, idest; as the time at which the mass of the 
main axial fragment has become smaller than the initial cloud mass by a factor of 1/e. The mass 
of the main axial fragment is evaluated as the mass of the gas that is physically connected to the 
main axial fragment by gas with p > 2pji. (Recall that the initial cloud is defined as the gas 
whose density is greater than 2piQ.) In practice, we identified the gas in the main axial fragment 
by drawing a density contour at a level of 2 pn and then integrating the density in the region 
surrounded by the contour line. 

In Tables 2 through 6, we summarize the cloud destruction times obtained from our numerical 
simulations. The cloud destruction time is longer for smoother clouds because both the Kelvin- 
Helmholtz and Rayleigh-Taylor instabilites grow more slowly in such clouds. For example, for 
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n = 2, tdest/^cc ~ 10 — 12, which is about three times larger than for the uniform case. On the 

other hand, tdcst/^cc depends only weakly on the initial density contrast, X; because the growth 
rates of the Kelvin-Helmholtz and Rayleigh- Taylor instabilities scale with tec (see the previous 
section). 

We note that our cloud destruction times of the uniform clouds tend to be smaller than those 
of KMC. For x = 10; '^e find i^est = 3.45 ice whereas they obtained i^est = 3.79 ice! for X = 100) 
we find ijest = 3.02 ice vs. their value of 3.90 tec- This is due to the fact that our cloud destruction 
time depends weakly on the threshold density of the main axial fragment, especially for large %. 
In this paper, the threshold density of the main axial fragment is set to twice the shocked ambient 
density, 2 pn, for all the models. For smaller values of the threshold density, the cloud destruction 
time increases: For example, at x = 100, tdcst — 3.6 tec for Pth = 1-5 Pii, and t^cst — 3.86 tec for 
Pth = 1-2 Pii- Thus, our cloud destruction times seem to approach those of KMC as the threshold 
density decreases, and we expect them to coincide with those of KMC if we were able to count in 
all the cloud material. 



7.2. Mixing Time: tmix 

Fluctuations with relatively short wavelengths, which have relatively large growth rates, are 
more effective at mixing cloud material with the the ambient gas than the long wavelength fluctu- 
ations that break up the cloud. To evaluate the efficiency of mixing, we define the mixing time, 
tmix, as the time at which the mass of the region whose density is greater than 2pji (md) reaches 
half the initial cloud mass (rricifi). 

In Figure 9, we depict the time evolution of the cloud mass normalized to its initial value, 
mci/rricifi- The cloud mass decreases more rapidly for smooth clouds than for uniform clouds at 
early times, whereas at late times, it decreases more slowly. This indicates that the lower density 
gas in the cloud envelope tends to be mixed more rapidly, but the high density gas in the central 
core tends to be mixed only after the formation of the slip surface, which is a prerequisite for the 
onset of the Kelvin-Helmholtz instability and which is retarded for a smooth clouds. For the n = 2 
models, the mixing time is determined by mixing of the envelope because of the large fractional 
mass of low density gas, whereas for models with larger n, the mixing time is determined by the 
KH instability. 

Because the mixing time depends on relatively small-scale fluctuations, one might expect it to 
depend sensitively on the spatial resolution of the computation. However, comparison with different 
spatial resolutions [R30, Reo, R120, and R240 for {n, x, M) = (8, 10, 10)] shows that the mixing time 
at i?i20 is within 8 % of that at R240, indicating that the mixing is reasonably accurate at -Ri2o- 
Also, the errors of the mixing time decreases nearly in proportion to the inverse of the cell number, 
which implies first order accuracy. 

Figure 10 shows how the cloud material is mixed with the ambient gas for two models: 
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(n, X, = (8,10,10) and (2, 10, 10). The histograms indicate the fraction of mass Am/mci,o 
within a corresponding density bin with a width of 0.2 /9co- For n = 8, about 50 % of the total 
cloud mass has density of ~ pcO the initial state. When the shock is transmitted into the cloud, 
the cloud density near the center increases by a factor of 4. At the same time, however, the shocked 
cloud starts lateral expansion, which reduces the density of the shocked cloud material. At t ~ tec, 
the distribution of the mass fraction is roughly flat over the range 0.4 pco < P < 4.4 pco- At this 
time, the cloud mass (the mass of gas with p > 2 pil) has dropped to about 85% of the initial cloud 
mass. Gas from the initial cloud that is at lower densities corresponds to envelope gas that has 
been shocked by both an incident shock and a reverse shock. Two shocks increase the entropy of 
the shocked gas above that that due to one shock. As a result the final density is only about twice 
the initial value (KMC), and some of the gas in the shocked envelope has a density below 2pii. 
Mixing due to hydrodynamic instabilities starts at t ~ 4fcc- As the cloud material mixes with the 
intercloud medium, the density of the mixed gas decreases with time, approaching the intercloud 
density with p ~ 0.4 pco- 

At t = 6.04 tec (Fig- 10 c), the mass distribution has two peaks, at p/pco ^ 0.7 and 2.5. The 

higher-density peak at p/pcO ^ 2.5 corresponds to the mass in the main axial fragment, which 
decreases with time. The lower-density peak at p/pcO ^ 0.7 corresponds to the mass that is mixed 
by the shocked ambient medium. The low density peak shifts towards lower density with time, 
gradually approaching the density of the shocked ambient medium, pn ^ 0.4 pcO- By t ~ 10 fee, 
most of the cloud material has densities lower than 0.8 pco- On the other hand, for n = 2, 
most of the mass (about 80 % of the total mass) is already contained in the lowest density bin, 
0.4 < Ap < O.GpcO- After passage of the shock, the density range of the cloud material widens 
immediately, although more than half the mass has a density less than p^o • The Kclvin-Helniholtz 
instability mixes cloud material with the ambient medium, thereby reducing the mass fraction of 
the high density material. By the time shown in panel (h) of Fig. 10, most of the cloud material 
has a density ~ 0.5 pco owing to efficient mixing. 

Our simulations demonstrate that the destruction and mixing of smooth clouds depend sen- 
sitively on the initial density distribution. In clouds with relatively steep density gradients, the 
core is almost completely destroyed by ^mix! for example, for n = 8 and % = 100, the core mass is 
0.10 niQ at t = tjnix (recall that in this case, the core mass is the mass with p > 0.51 Pco)- Thus, 
at this point the high density core is almost destroyed, and cloud mixing proceeds very efficiently 
throughout the entire cloud. By contrast, significantly more of the core survives for clouds with 
shallow density gradients: mco(imix) = 0.29mco,o for the model with n = 4 and x = 100, and 
mco{tmix) = 0.25mco,o for the model with n = 2 and x = 100. For the case of a shallow density 
gradient, the cloud envelope is stripped from the main body of the cloud by t = tmix, but much 
of the core remains intact. In such a core, the velocity dispersion is also relatively small, which 
indicates that random motions produced by the shock-cloud interaction arc directly associated with 
cloud destruction, and do not significantly affect the cloud core before the core is destroyed. 
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8. Vorticity 

The interaction of a shock wave with a cloud generates powerful vortex rings, which play an 
important role in the destruction of the cloud. For uniform clouds, KMC constructed a simple 
analytic model for vorticity production by the shock-cloud interaction that was in good agreement 
with their numerical results. In this section, we extend their model to the smooth clouds, and 
compare it with the numerical results presented in the previous sections. 



dt I i 



8.1. Vorticity Production: Analytic Model 

In the course of the shock-cloud interaction, vorticity is produced by the baroclinic term (the 
right-hand side of eq. [28]) in the vorticity equation for an inviscid fluid, 

— + Vx{u}Xv) = VPxV\-\. (28) 

Integrating the above equation over a cross section A, we obtain the increase rate of circulation 

VP X V^J -dA . (29) 

According to KMC, the vorticity production can be classified into four components. Two are 
associated with the cloud-intercloud boundary and are produced by the initial passage of the shock 

(Fshock) and by the subsequent post-shock flow (Tpost)- The third is associated with the triple 
points associated with the Mach reflected shocks behind the cloud (Fring)- The fourth is the 
vorticity produced in the cloud (Fdoud); which is smaller than those in the other components by a 
factor of order x~^/^. As shown below, the vorticity production in smooth clouds is qualitatively 
similar to that in uniform clouds. 

First, we consider the circulation produced by the initial passage of the shock, Fshockj defined 
as the circulation generated by the time at which the shock sweeps over the cloud and converges 
on the z-axis. For the uniform case, this component is generated only in a thin layer between the 
cloud and the intercloud medium. In contrast, for smooth clouds, the layer where this component 
is produced is extended in the radial direction because of the smooth cloud boundary. Indeed, as 
we shall see below, for smooth clouds such as those with n = 2 and x = 10; this circulation is 
produced throughout most of the cloud envelope, not just at the cloud surface. The time for the 
shock to sweep over the cloud is tsweep ^ Scto/fft. For uniform clouds, this estimate of the time scale 
coincides with that of KMC. KMC estimated the increase rate of Fghock as 

dF 3 /vhCos6\ / dz/ cos 6\ , i/.,, So, i/n, , , 

where 9 is the angle at which the shock impacts the surface and h is the thickness of the shear 
layer, which is approximately the grid-cell size for uniform clouds. A strong shock is also assumed 
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(see eq. 7.8 of KMC for more details). KMC allowed for a reduction in the shear due to the motion 
of the shocked cloud by inserting a factor (1 — x"^''^) iii their final result for rghocki here we have 
included this factor in the original equation as well. 

Equation (30) indicates that the rate of increase of Tshock is independent of the thickness of 
the shear layer. Thus, the fact that smooth clouds have a shear layer with a thickness h /ibdj the 
scale height in the boundary (eq. 5), should not affect the growth of the circulation. Thus, we find 

Tshock - ^ ^ *sweep - ~Vbao{l - X"^^^) ■ (31) 

In the next subsection, we show that this simple model reproducesour numerical results reasonably 
well. 

Next, we consider the vorticity production due to the interaction of the shocked cloud with 
the post-shock flow, Tp^st- When the initial shock passes through the cloud and the slip surface is 
formed at the upstream side of the cloud, the vorticity continues to be generated at the slip surface 
by the baroclinic term, which is non-zero when the density and pressure gradients arc not parallel, 
and the equation of state is not barotropic. The pressure is approximately constant across the slip 
surface, but it varies along the slip surface, being a maximum at the stagnation point at the front 
of the cloud and smaller along the sides of the cloud. On the other hand, density varies primarily 
across the slip surface. Since the slip surface resembles the surface of the uniform cloud, we adopt 
KMC's equation for the growth of the post-shock circulation: 

APA {-] ^ ^ — v'^ . (32) 



(it \pj pii 2 " 

Then, the circulation produced by the post-shock flow can be estimated as 

Tpost - -^^Stdrag - (^^^) (1 - X-'^rv,r^o , (33) 

where we have approximated the cloud velocity as v'^. = 6xp(— t/tdrag) with v'^ = [2vi,/{'y + 
!)](! — x~^^^)- Here we have followed KMC in adopting a simpler approximation for the cloud 
velocity than that given in §2.3 (eqs. [24] and [25]). 

Third, we consider the vorticity production associated with the triple points that are formed 
by the interaction of the initial shock and the reflected shock behind the cloud (Tring)- Following 
KMC, we estimate the circulation of the supersonic vortex ring as 

3 , X 

Tring = -^Vbtto . (34) 

Note that this component has a positive sign, which is the opposite of Fghock and Tp^st (see §7.1 of 
KMC for more detail). 

Finally, the vorticity produced in the core Fcorei is expected to be smaller than Fghock and Fpost 
by a factor of order x~^/^. Since this is small, we follow KMC in not attempting to model this 
component. 
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Recently, Kornreich & Scalo (2000) constructed an analytic formula for the vorticity generated 
by the shock-cloud interaction for smooth clouds, which corresponds to our Tshock- They assumed 
that (1) the post-shock flow is steady and parallel to the z-axis, and (2) the velocity induced by 
the shock is inversely proportional to the square root of the pre-shock density. In order to make a 
comparison with our analytic model, we derive the circulation estimated from their model. Prom 
their eq. (40), 



Figure 11, shows the time evolution of the total circulation, F, for several different models. For 
all the models, the evolution is qualitatively similar: F increases monotonically until the intercloud 
shock reaches the rear of the cloud, then decreases with time as the shock transmitted into the 
cloud interacts with the reflected shock on the z-axis. Thereafter, the total circulation increases 
again. The first peak corresponds to Fghock because at this time the incident shock has passed 
over the cloud and reached the z-axis at the rear of the cloud. The subsequent decrease in F 
corresponds to the positive contribution due to the supersonic vortex ring, Fring. The interaction 
of the shocked cloud with the post-shock flow continues to generate the circulation Fpost, so that 
the total circulation subsequently increases with time. 

For clouds with relatively large density gradients [ (n, x) = (8,10) and (8,100)], the values 
of Fghock/'^co^b are predicted from eq. (31) as —1.77 and —2.58, respectively, which are within 
about 15% and 30% of the peak values of the total circulation of —2.01 and —3.78, respectively. 
Kornreich &: Scalo (2000)'s model (see eq. 35) predicts Fks = —1.30 and —2.18 for the two cases, 
which is somewhat less accurate than our model. Our numerical results show that the circulation 
associated with the supersonic vortex ring is about 1.02 for (n, %) = (8, 10), which is close to our 
analytic estimate of 0.866. For (n, %) = (8,100), the numerical and analytic results are 1.40 and 
0.953, respectively, again in good agreement. The analytic model predicts that the circulation due 
to the post-shock flow is Fpost — — (0.48, 4.4)rcoWft for the models with x = 10 and 100, respectively, 
which is in reasonable agreement with our final results. At the end of the computation, the total 
circulation is continuing to increase. However, we expect the rate of increase to decrease with time 
as the cloud becomes comoving with the post-shock flow so that the circulation will not increase 
much thereafter. Lower resolution calculations confirm this prediction. 

The components of Fpost and Fring can be clearly seen in the spatial distribution of the circula- 
tion shown in Fig. 11, in which the vorticity is integrated from the bottom of the computation box 




(35) 



We compare this to the numerical results below. 



8.2. Comparison With Numerical Simulations 
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up to The cumulative circulation at t = 2.86 shown in Fig. 12f, decreases with z, reaching a 
minimum value of —2.3 VcoVb, which corresponds to Tshock- As z increases further, the circulation 
increases due to the contribution of the supersonic vortex ringj Tring — I'O^co^fe* For clouds with 
relatively steep density gradients (like n = 8), the spatial distribution of vorticity generated by 
the shock-cloud interaction is qualitatively similar to that of the uniform cloud (see Fig. 16 of 
KMC), i.e., immediately after the initial shock passage, the vorticity generated by the interaction 
is confined to a thin layer between the cloud and intercloud medium (Fig. 12b). In the subsequent 
evolution, the vorticity is concentrated in the slip surface, as can be seen in Fig. 12c. 

For clouds with small density gradients, the analytic model continues to work well for the 
circulation produced by the shock. In this case, to estimate the circulation from the analytic 
model, we adopt the reduced initial cloud radii used in §5.2. For {n,x) = (2, 10), our simulations 
show that the first peak in the circulation is — S.lOrcofb, compared with Fghock = — 3.41rcoi'6 from 
our analytic model. For (n, x) = (2, 100) the corresponding values of the circulation are —12.6rcoVb 
and —lA.lrcoVb, respectively. Again, the analytic predictions are within about 12 % of our numerical 
results. In contrast, equation (34) overestimates the vorticity production by the supersonic vortex 
ring for the model with n = 2 and x = 10 by almost a factor 3 (our numerical results give 
Tring — O.SOlrcoi'b, whereas the analytic model predicts Fring = 1.66rcof6). In this model, the 
supersonic vortex ring is not as powerful as for the models with larger density gradients. As 
mentioned in §4.2, the vortex ring propagates more slowly than the triple point that moves with 
almost the same velocity as the incident shock Vb- This slower propagation speed of the supersonic 
vortex ring appears to be special to the n=2 and x = 10 model. For example, for the model with 
n = 2 and x = 100, a powerful supersonic vortex ring is formed after shock passage. The numerical 
result shows Fj-ing — 5.50rcoVb, whereas the analytic solution predicts Frmg = 5.25rco^^6, about 5 % 
less than our numerical result. 

For smooth clouds with n = 2, the analytic model of Kornrcich & Scalo (2000) predicts 
^Ks/'/'coVb = —3.98 and —14.5 for x = 10 and 100, respectively. These values are in reasonable 
agreement with the numerical results given above, although once again the error is somewhat larger 
than for our analytic model. 

For smooth clouds, the vorticity generated by the initial interaction with a shock extends 
throughout the cloud envelope (see Fig. 13b). However, this radially-extended vorticity shifts 
towards the slip surface with time. Eventually, the vorticity is concentrated in the slip surface, 
and is subsequently converted to random motions of small fragments as small-scale fragmentation 
proceeds. 

We note that the evolution of the total circulation for the uniform clouds (dotted curves in 
Fig. 11) are apparently different from those of KMC (Figs. 13 and 15 of KMC) after the incident 
shock reaches the z-axis at the rear of the cloud. This discrepancy occurs because the contribution 
of the supersonic vortex ring is always included in our calculations, whereas KMC did not include 
it after the ring left the high-resolution part of the grid. If the contribution due to the supersonic 
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vortex ring is removed from our results, the evolution of the circulation is in good agreement with 
KMC's results. 

In summary, the analytic model constructed in §8 agrees well with our numerical results. 
Recently, Kornreich & Scale (2000) claimed that the vorticity production in smooth clouds is not 
related to the baroclinic term, but is purely kinematic in origin; the preshock density distribution, 
which varies radially through the cloud, naturally generates the velocity gradient in the postshock 
cloud. However, the agreement between our theory and our simulations have clearly demonstrated 
that the vorticity is indeed produced by the baroclinic term in eq. (28) during the shock-cloud 
interaction. Their analytic model for vorticity production by the shock is somewhat less accurate 
than ours, since they did not consider the vorticity produced in the post-shock flow or in the 
supersonic vortex ring. 

9. Discussion 

9.1. Effect of a Smooth Cloud Boundary on Cloud Morphology: Implications for 

Observations 

In this paper, we have studied the interaction of shocks with interstellar clouds with smooth 
density profiles. We showed that after several cloud crushing times, a slip surface with a significant 
density jump develops in front of the core of the shocked cloud, irrespective of the initial density 
distribution. Thus, the cloud is deformed into a head-tail structure before subsequent cloud de- 
struction (see e.g.. Figs. 2c and 3d). This suggests that the density profiles of observed shocked 
clouds should have significant density jumps in the direction of the source of the large scale shock. 
One might think that the effects of the initial density profile on cloud morphology would be ob- 
scured as the slip surface develops, but in fact information on the density profile remains imprinted 
on the shocked cloud until significant instabilities develop. 

Figure 14 shows the time evolution of the ratio of the masses mo. 5 and mo.i for the models 
with (n,x,M) = (2, 100, 10) and (8, 100, 10), where mo. 5 and mo.i represent the masses contained 
within the regions with p > 0.5 /9max and p > O.lpmax, respectively. Prior to cloud destruction, 
the mass ratio mo.s/mo.i remains approximately constant, but varies significantly with n. For 
example, in the early stages of the shock-cloud interaction, the mass ratio stays at mo.i/mo.5 > 0.5 
for clouds with relatively steep density gradients (n > 8), whereas it is only about 0.1 for smooth 
clouds (n 2). This suggests that observations of this mass ratio may be able to reveal the initial 
density profiles of the preshock clouds provided the clouds have not entered the shock destruction 
stage; the initial mass ratios 71^^ listed in Table 1 could be directly compared with the observed 
values. 

Our simulations can be compared with the observations of young and middle-aged SNRs such 
as Cas-A and Cygnus Loop. Two important caveats in making any such comparison are that (1) 
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our simulations are non-magnetic, whereas magnetic fields are believed to be important in the 
interstellar medium, and (2) almost all of our simulations are adiabatic, whereas radiative losses 
are important in shocked interstellar gas. Recently, Patnaudc et al. (2002) and Danforth et al. 
(2001) performed detailed investigation on the shock-cloud interactions observed in the Cygnus 
Loop (the XA and southwestern region of the Cygnus Loop). Their observations do not show any 
signs of the Kelvin-Helmholtz instability at the observed cloud surface, and thus they concluded 
that the interaction of the shock and cloud is in the initial transient stage. However, if the prcshock 
cloud had a smooth density gradient at the cloud boundary, then the Kelvin-Helmholtz instability 
does not grow in the shock compression or the early post-shock expansion stages. If the mass ratio 
mo.s/mo.i (or an equivalent one) could be estimated from the observations, then one might be able 
to infer something about the structure of the initial cloud. 



9.2. Velocity Dispersions of Shocked Clouds 

As shown above, the shock-cloud interaction produces vorticity inside the cloud. Much of the 
vorticity is concentrated at the slip surface, which is subject to the Kelvin-Helmholtz instability, 

and as a result the shocked cloud is distorted, producing vortex rings. Vortical motions generate 
substantial velocity dispersions. Recently, velocity dispersions generated by the shock-cloud inter- 
actions have attracted a great deal of attention as a mechanism for replenishing turbulent motions 
in the ISM (e.g., Kornreich & Scalo 2000; MacLow & Klessen 2003). Here, we discuss how velocity 
dispersions generated by shock passage evolve with time. First, we compute the synthesized spectra 
of turbulent motions integrated along the z-axis (see e.g., Falgarone et al. 1994), and then calcuate 
the dependence of the velocity dispersion on region size, the so-called linewidth-size relation. 



9.2.1. Synthesized Spectra of Shocked Clouds 

Figures 15a through 15c show the evolution of the synthesized spectra of a shocked cloud 
with {n,x,M) = (8,10,10). The synthesized spectra have been computed by averaging density- 
weighted histograms of the line-of-sight velocity over the entire cloud component. As mentioned 
in the next subsection, the velocity dispersions are typically subsonic for adiabatic shocks, and 
thus thermal broadening would smear the motions of small cloud fragments. Thus, to understand 
how the shocked cloud fragments due to hydrodynamic instabilities, we focus on the non-thermal 
motions. The line-of-sight is assumed to be parallel to the z-axis. At t = 4 tec, the synthesized 
spectrum has a fairly smooth profile with a peak at O.dvn, where vn is the velocity of the shocked 
intercloud medium, and a tail at higher velocities. At this stage, the cloud has not been destroyed, 
although it has been deformed into a head-tail structure (see Fig. 2d). The velocity component near 
the peak corresponds to the main body of the shocked cloud, while the higher velocity component 
corresponds to the tail of the shocked cloud. By t = 8 tcc„ the cloud has broken into small fragments, 
mainly due to the Kelvin-Helmholtz instability (see Fig. 2f). Each fragment has a different velocity. 
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which produces many spikes in the synthesized profile. The velocity dispersion of these fragments is 
estimated to be of the order of 0.1 Vf,, somewhat less than the sound speed behind the cloud shock 
(~ O.lSvb). The synthesized spectrum has a profile similar to panel (b) for 6tcc ^ ^ ^ lOtcc- By 
t = 13.2 tec (Fig- 15 c), these spikes have become fainter and the peak velocity is shifting towards 
the post-shock ambient velocity as the cloud is accelerated and mixed with ambient medium by the 
interaction with the postshock flow. The component near the peak velocity has a velocity dispersion 
^ O.Q5vb. The synthesized spectrum also has a plateau extending towards small velocities, which 
corresponds to the cloud components near the z-axis that were destroyed recently. This plateau 
becomes smaller with time. 

Figures 15d through 15f show the evolution of the synthesized spectra for a shocked cloud with 
{n,x,M) = (2, 10, 10). The overall evolution of these synthesized spectra is similar to that for the 
n = 8 case. Prior to the onset of the Kelvin-Helmholtz instability, the spectrum is more-or-less 
continuous and has a peak at low velocities, extending towards higher velocities. Bjt = 12 tec j the 
cloud has begun to be destroyed and a number of spikes appear in the spectrum. The synthesized 
spectrum has a profile similar to panel (e) in the time interval of 10 tec ^ ^ ^ 20tcc- At the late 
stage of cloud destruction, the peak of the spectrum shifts towards the velocity of the post-shock 
ambient medium. The evolution of the spectrum is slower than that of the n = 8 case because of 
the smaller relative velocity between the cloud and intercloud medium. 

The dissipation of random motions seems to be slow compared to the cloud-crushing time. 
The velocity dispersion decreases by about 30 % from the maximum after t = tdest for n = 8, 
while for n = 2, the velocity dispersion does not decrease significantly after t = idest- On the other 
hand, for the uniform cloud (KMC), the velocity dispersion decreases by about 40 % by the end of 
computation (10 tec)- This reflects the fact that smoother clouds tend to evolve more slowly. The 
velocity dispersions would decrease more if further evolution is followed for smoother clouds. In 
the next subsection, we discuss how dissipation of random motions depend on the size, based on 
three-dimensional simulations. 

9.2.2. The Linewidth-Size Relation (3D simulation) 

Observations of the ISM indicate that interstellar clouds have nonthermal internal motions that 
are inferred to be supersonic turbulence. The observed velocity dispersions are found to scale with 
region size i? as a power law, Av ~ R^, with p = 0.4 — 0.6 (e.g., Larson 1981; Myers 1987; Dame et 
al. 1986), although the power index p can depend on the environment (e.g., Caselli & Myers 1995). 
In this paper, we have shown that the interaction of shocks with interstellar clouds can generate 
random motions in the cloud material, but only by completely destroying the clouds. Here, we 
calculate the dependence of the velocity dispersion on region size, the so called linewidth-size 
relation, from our numerical results, and discuss whether random motions generated by the shock- 
cloud interaction agree with the observed scaling law. Because the characteristics of hydrodynamic 
turbulence depend strongly on the spatial dimension, we perform a three-dimensional calculation 
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with an intermediate spatial resolution Rqq for the nonradiative cloud with (n, %, M) = (8, 10, 10). 
We confirmed that all the global quantities in the three-dimensional calculation coincide with those 
of the two-dimensional calculation with the same spatial resolution to within 10% in the time 
interval < t < lOicc- 
In Fig. 16, we represent the linewidth-size relation for the model with (n, M, x, 7) = (8, 10, 10, 5/3) 
at three different times: t = 6 tec, Sico and 10 ice- The abscissa denotes the beam size, while the 
ordinate denotes the velocity dispersions of the cloud material. The figure shows the results of 
observing the shocked cloud from two different directions, (1) along the z-axis (perpendicular to 
the incident shock plane; shown by solid curves), and (2) perpendicular to the z-axis (in the plane 
of the incident shock; shown by dotted and dashed curves). In the latter case, we consider two dif- 
ferent positions for the center of the observing beam: at the center of gravity, zq, and downstream 
from that point, zq + Tco, where the cloud is more disrupted than at z = zq- 

By t = 6 tec 5 the cloud has begun to experience significant destruction due to hydrodynamic 
instabilities, but at this stage the main axial core is still recognized on the z-axis (sec Fig. 2). 
The solid and dotted curves have peaks at ~ 0.3rco and ~ 0.6rco, respectively. These peaks 
correspond to the components near the main axial core. In the downstream observing beam, much 
of the cloud is mixed with the intercloud material, and the velocity dispersion (dashed curve) 
is almost independent of scale. This part of the cloud has been destroyed due to the nonlinear 
growth of small-scale Kelvin-Helmholtz instabilities. This indicates that the turbulent motions are 
generated mainly in the cloud envelope that has been destroyed, and in the cloud core, the turbulent 
generation delays. By t = 8tcc; the velocity dispersions in all the observing beams become relatively 
fiat, indicating that the entire cloud is being destroyed due to the Kelvin-Helmholtz instabilities. 
Thereafter, the velocity dispersions in the small scale dissipate, as can be seen in Fig. 16(c). As 
a result of dissipation, at the time of t = 10 tec; if the velocity dispersions are approximated by 
power-laws as Av ~ Rf, the power index is estimated as p = 0.2 — 0.3, although the power index 
depends on the direction and beam position. At later times, p increases somewhat since the small- 
scale fluctuations damp out more rapidly. Many numerical studies show that interstellar turbulent 
motions dissipate in a dynamical time, L/Sv, where L is the size of the region (Elmegreen & Scalo 
2004). At t ~ 8 tec, we find L/Sv ~ 0.53tcc and 2.9icc at L/rco ~ 0.1 and 1.0, respectively, where 
we have averaged over three different beams. The velocity dispersion at small scales, r ~ O.lrco, 
dissipates significantly from t = 8tcc (Fig- 16b) to t = 10 tec (Fig. 16c), whereas at large scales, 
r « Irco, turbulent dissipation is less efficient. It appears that there has not been enough time to 
establish a self-consistent turbulent cascade to maintain the turbulent motions on small scales. 

Our simulations show that at the early stage of cloud destruction, small-scale fluctuations 
dominate, which makes the linewidth-size relation flat. Subsequently, these small-scale fluctuations 
tend to damp out in a few cloud crushing times, leading to a linewidth that increases with scale. 
The linewidth-size relations obtained from our simulations are time-dependent and have not reached 
their asymptotic values. At late times, the cloud material dissolves into the ambient medium due 
to mixing, and therefore, it is difficult to distinguish them well. The turbulent cascade could also 
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influence the evolution of the hnewidth-size relations, although our resolution is inadequate to 
study this precisely. Sytine et al. (2000) performed convergence tests on decaying turbulence with 
periodic boundary conditions and suggested that at least 1024^ grids are needed to capture the 
characteristics of turbulent fields accurately. However, in our 3D simulation with Rqq, the effective 
radius of the shocked cloud expands by a factor of two by i = tacst, so that the shocked cloud 
is covered by only about 256^ grids at the epoch of cloud destruction. Thus, to discuss the time 
evolution of the linewidth-size relation more accurately, higher resolution runs may be needed. 

9.3. Implications for Interstellar Turbulence 

The ISM is observed to be highly turbulent. Recently, Hcilcs Sz Troland (2003) found that 
the mass- weighted velocity dispersion is 7.1 km s^^ for cold neutral medium (CNM) clouds and 
11.4 km s^^ for warm neutral medium (WNM) clouds. The one-dimensional rms Mach number of 
internal motions for individual CNM clouds is ~ 1.7. The mean temperatures of the CNM and 
WNM are measured to be ~ 70 K and ~ 8000 K, respectively. Under the assumption of pressure 
equilibrium between the CNM and the WNM, the typical density contrast between CNM and WNM 
clouds is of order x = 100. Heiles & Troland (2003) argue that the CNM clouds are organized into 
large sheets. More recently, Stanimirovic & Heiles (2005) have found very small CNM clouds, with 
column densities ~ 10^^ cm~^. In this subsection, we discuss how our results on the shock-cloud 
interaction might be related to these observations. 

Interstellar shocks arc produced by several events such as supernova explosions and stellar 
winds from massive stars. Among them, the injection rate of energy into the ISM through super- 
novae is expected to exceed those of other sources (e.g.. Castor 1993). Thus, we concentrate on 
the effects of supernova remnants (SNRs) and attempt to estimate the magnitude of the generated 
velocity dispersions. Most of the shock-cloud interactions arc expected to be at low Mach numbers, 
because the shock speeds of SNRs rapidly decrease with expansion, although the typical Mach 
number is highly uncertain. In the following, we adopt the typical Mach number of M ~ 3 as 
a representative value. For shocks at M ~ 3, the distance from the origin of a blast wave Rc is 
estimated to be of the order of 100 pc (Cioffi, McKee &; Bertschinger 1988). In this case, the small 
cloud approximation is applied for clouds with radii of the order of 1 pc or less, (see §2.2). This 
size is comparable to that of typical CNM clouds observed by Heiles & Troland (2003). 

As shown above, the velocity dispersion induced by the shock-cloud interaction is of the order 

of 

6v ~ 0.074vb = 2.2M3X2^Cco,5 kms"^ , (36) 

for radiative clouds with 7 = 1.1, where M3 = M/3, X2 = and Cco,5 = Cco/10^ cm s""^. The 

above one-dimensional velocity dispersion [Sv'^ = {2Sv'^ + 5v'^)/3] is estimated as the mean value 
of the n = 8 and n = 2 clouds at t = tm- For T = 70 K, the one-dimensional velocity dispersion 
is computed as about 1.4 km s~^. This magnitude of turbulence is comparable to the observed 
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internal velocity dispersions of CNM clouds, Sv ~ 1.1 km s ^, (Heiles & Troland 2003). 

When one observes the ISM, it is difficult to directly distinguish the structures along the line- 
of-sight from observations. Approximating the thickness of the CNM to be 6 = N^ji/ubi, Heiles & 
Troland (2003) infer that the CNM is in thin sheets with typical aspect ratios in the range of 100 
— 300, where A^hi and nni are the column density along the line-of-sight and the volume density, 
respectively. If one observes a shocked cloud from the side and models it as a sheet of length 2c, 
width 2a, and thickness b, then the thickness b is evaluated as 

riHi (2a X 2c) (pi) 3oc ' ^ ^ 

where we consider isothermal shocks. Then, the inferred aspect ratio is given by 

2c/6 = 6ac^M'^/{nal) . (38) 

For a shocked cloud, the above aspect ratio tends to be large and increases with time because the 
length c increases with time. For example, for a shocked cloud with n = 8, x = 100, and M = 3, 
the aspect ratio 2c/6 reaches about 2000 at t ~ tdcst- This implies that some of the observed 
CNM clouds with large aspect ratios could be shocked clouds torn up into small fragments, instead 
of thin sheets. Small cold clouds recently observed by Stanimirovic &; Heiles (2005) could be 
related to this case. Since we don't take into account the magnetic fields, we cannot compare our 
result quantitatively with observations. The omission of the magnetic fields most likely leads to 
overestimation of the shocked density and therefore the aspect ratio. Although shocks in the real 
ISM are not isothermal as we assumed above, deviations from isothermality behind the shocks 
would cause only quantitative, not qualitative, changes in the result, provided that the cloud was 
initially CNM. 

In our simulations, the vorticity generated by the shock-cloud interaction is converted into 
random motions of small cloud fragments produced by the Kelvin-Helmholtz and Rayleigh- Taylor 
instabilities. In that case, the CNM may consist of small shreds in a coherent flow. We note 
that this mechanism of turbulence generation is different from that of Kornreich &; Scalo (2000), 
who proposed that the vorticity produced by the shock-cloud interaction is confined inside the 
shocked cloud and eventually converted into the turbulent motions observed in the interstellar 
clouds. They assumed that the smooth clouds are not destroyed significantly by shocks and that 
the density profile of the cloud is more or less maintained after shock passage. They suggest that 
the velocity dispersions generated by shock-cloud interaction can be scaled as Av/vb = n^/^(pi — 
po)/{2pi)x~^'^ oc x'"^'^. However, our simulations have demonstrated that the density profile in 
the cloud changes significantly with time, and within a few cloud-crushing times a slip surface 
forms in front of the core of the cloud. In fact, the velocity dispersions are estimated as O.lvh, 
which is more or less independent of Furthermore, prior to complete cloud destruction, turbulent 
motions produced by the shock-cloud interaction are not induced in the cloud core, but in the cloud 
envelope. For example, for the model with {n, x, M) = (2, 100, 10), the circulation in the cloud core 
(Fco), which has a density greater than 1.96 Pco [~ 4:po{r = Tco)] is ^co/fcoVb = 0.089 at t ~ 1.5 tec- 
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This value is only about 0.7 % of the total circulation {Ttot /fcoVb ^ 12), indicating that most of 
the vorticity is in the cloud envelope. (If we use eq. 40 of Kornreich & Scalo (2000), which gives 
the distribution of vorticity as a function of radius, the circulation generated by the shock passage 
is estimated as Tks / '"co Vb = 0.039 in the cloud core of r < rco and TKs/rcoVb = 14.5 in the whole 
region of r < oo. These estimates from the KS theory also show that the shock-cloud interaction 
is ineffective at inducing turbulent motions in the cloud core.) For other models with larger n, the 
ratios of Fco to Ftot are much smaller than for n = 2. 

In summary, the shock-cloud interaction can produce the magnitude of the velocity dispersions 
observed in the typical CNM clouds with the size of ~ 1 pc. However, turbulent motions produced 
by the shock-cloud interaction are directly associated with cloud destruction and, in the absence of 
a magnetic field at least, do not significantly affect the cloud core, contrary to Kornreich & Scalo 
(2000) 's view. The interaction of a shock with a cold H I cloud is expected to lead to the production 
of a spray of small H I shreds which have a coherent velocity distribution. Such structures could 
be related to the small cold clouds recently observed by Stanimirovic &: Heiles (2005). However, 
we cannot make a realistic estimate for the size of the smallest shreds to compare with their 
observations, because our calculations do not include magnetic fields and thermal conduction (see 
e.g., Koyama &; Inutsuka 2002 for the effect of thermal conduction). 



9.4. Implications for Star Formation 

Under some circumstances, compression due to interstellar shocks is expected to make a 
shocked cloud self-gravitating even if self-gravity is unimportant in the preshock cloud, thereby 
triggering star formation (Elmegreen Sz Lada 1977). Foster &: Boss (1996) and Vanhala &; Cameron 
(1998) numerically studied the interaction of radiative shocks with self-gravitating spherical clouds, 
showing that gravitational collapse can be induced by shock compression when radiative cooling 
is efficient in the shocked clouds. In particular, Foster &: Boss (1996) studied the interaction of a 
critical Bonnor-Ebert sphere and a thin pulse-like wave that mimics an interstellar shock or wind, 
under the assumption that the gas is isothermal. They found a critical momentum of this wave 
above which gravitational collapse is induced after the interaction. Unfortunately, their result can- 
not be directly compared with ours because their initial conditions are different from ours; e.g., 
their initial wave is not satisfied with the shock jump condition. Vanhala & Cameron (1998) per- 
formed SPH simulations of the interaction of a planar shock and a spherical, centrally concentrated 
molecular cloud core, taking into account the effect of radiative cooling. They found that shocks 
with intermediate speeds of 20 — 45 km s^^ are capable of triggering gravitational collapse since 
the effective 7 of the shocked cloud remains close to unity. In contrast, stronger shocks inhibit grav- 
itational collapse because the main coolants such as CO and H2 are destroyed for velocities above 
about 40-50 km s^-^, driving the effective 7 above 4/3. This critical velocity for shock destruction 
of the molecules is consistent with that of Draine & McKee (1993). More recently, Hcnncbclle et 
al. (2003, 2004) studied the evolution of molecular cloud cores compressed by external pressure and 
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demonstrated that a compression wave triggers collapse from the outside in. 

Here we first show that only radiative shocks can trigger gravitational instability in shocked 
clouds. We show that turbulence behind the shock has a strong stabilizing effect. Shocked clouds 

are compressed into a sheet-like configuration, and we analyze the condition for gravitational in- 
stability for this case following the treatment of Elmegrccn & Elmegreen (1978). Since we have not 
included the effects of magnetic fields in our simulations, we neglect them here. Magnetic fields 
tend to stabilize against gravitational collapse, so the results below give necessary but not sufficient 
conditions for gravitational instability. 



9.4- 1- Maximum critical mass 
Define the dimensionless cloud mass as (McKee & Holliman 1999) 



H = TTT-cl 



(39) 



where pa and Ccs arc the density and adiabatic sound speed at the surface of the cloud. A stable 
isothermal cloud has p < 1.18, the value for a Bonnor-Ebert sphere. More generally, a stable cloud 
has fi < ficr, where is a number that, for polytropes (P oc p'^^), depends on both the polytropic 
index 7^ and the adiabatic index 7 (McKee & Holliman 1999). The maximum mass of a cloud in 
hydrostatic equilibrium (i.e., the cloud mass corresponding to jUcr in eq. 39) is 

/p _4\3/2 

rricT oc f -J • pi -M , (40) 

The first factor is proportional to the exponential of the entropy, and therefore always increases 
in shocks unless they are sub- isothermal (i.e., the post-shock temperature is less than the pre- 
shock temperature). Since the density also increases in shocks, we conclude that shocks increase 
the critical mass — i.e., they stabilize clouds against gravitational collapse — for 7 > |. This is not 
surprising since clouds that are approximately locally adiabatic are stable against collapse for 7 > | 
in any case. 

For 7 < I , it is possible that the effect of the increased density might overcome the effect of the 
increased entropy in the expression for the critical mass. In fact, in the real ISM, radiative cooling 
can sometimes reduce the entropy of the postshock gas below the prcshock value (a sub-isothermal 
shock). In that case, the above critical mass is likely to decrease below the preshock value, if 
turbulent motions generated by shock passage do not contribute to cloud support significantly. 

In order to treat sub-isothermal shocks, we write the shock jump conditions as 

pivi = povo , (41) 
pM + eCl,^o) = Po{v^o + CL,o) , (42) 
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where Qso = {P/ pY^"^ is the isothermal sound speed and e = C^^ i/C'?^ q. Thus, e = 1 corresponds 

to an isothermal shock and e < 1 to a sub- isothermal shock; normal adiabatic shocks have e > 1. In 
shocked CNM clouds, the temperature might decrease to about 10 K if molecules can form behind 
the shock; in that case, e ~ 0.1. ^From the above relations, the post-shock density is 

M? 

Pi = -^Fpo , (43) 
where Miso = vq/ Ciso, o is the isothermal Mach number and 

F= V ■ _ If^ _ (44) 



The factor F approaches unity either for e ^ 1 or for Miso ^ oo and e/M^^ — > 0. In general, 

the factor F is about unity for |1 — el/M^^ <C 1, which is true for sub- isothermal shocks provided 
the isothermal Mach number Miso is not small; e.g., for e = 0.1 and Miso = 3, F ~ 1.1. Thus, 
the post-shock density can be approximated by pi ~ M^^^po/e for isothermal and sub-isothermal 
shocks. 

We can obtain an upper bound on e such that the shocked cloud is gravitationally unstable 
by assuming that the shocked cloud material is as localized as possible-i.e., it is spherical. Recall 
that /X is the ratio of the cloud mass to the critical mass (eq. 39), so that 



P_ 
Po 




(45) 



where /io is the value of /x for the pre-shock cloud. Now, Pi ~ M^^^Pq ^'iid, for P ~ 1, pi ~ 
(M^q/€)po from equation (43). A necessary condition for gravitational instability is p, > pci, which 
requires 

e < (MoMiso/Mcr)^/^. (46) 

Under the assumption that the pre-shock cloud is gravitationally stable, we have po < pcr, so the 
value of e in this equation satisfies the condition that P ~ 1 (i.e., e/M?.^ <g; 1) for Miso > 3. This 
equation suggests that arbitrarily strong radiative shocks can induce gravitational collapse, but two 
effects prevent this: non-thermal motions behind the shock (considered immediately below) and 
inadequate time for the gravitational instability to develop (considered in §9.4.2. 

The shock-cloud interaction generates nonthermal motions behind the shock, both random 
and systematic, that tend to stabilize the cloud against gravitational collapse. In a radiative shock, 
these motions will dominate the thermal motions behind the shock so that equation (45) can be 
rewritten as 

/ r< \ 3 / \ 1/2 

p _ ( Oc,iso,o W _ (47) 



Mo V "^^i / VPo 
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Since Svi increases with scale, it must be evaluated on the same scale as /x/^o- Our calculations are 
two-dimensional, so we estimate the one-dimensional velocity dispersion as 5v'^ = i25v'^ + 5v1)/2>. 

1 /2 

Behind the shock, we express 5v as a fraction of the shock velocity, 5v\ = KVb ~ k (x) ' vq, so that 

^ = """" (48) 

where we used equation (43) to evaluate the density jump. As remarked above, 6vi and therefore k 

increase with scale. When averaged over the entire shocked cloud, our numerical calculations give 
AC ~ 0.04 — 0.07 for t ~ (1.5 — 2)tcc for a uniform cloud with x = 100. Therefore, in typical CNM 
clouds, which have x ~ 100, the condition for instability {fj, > ficr) restricts the Mach number to 
be 

r-l/l 1/2 

^iso < ^ f '' , (49) 

where k_i = k/0.1 and X2 = x/100. For isothermal and sub-thermal shocks (0 < e < 1), we have 
1 < F < 2. We conclude that turbulence behind the shock stabilizes gravitational instability on 

the scale of the cloud for dense clouds {x2 ^ 1) unless the shock is relatively weak. Even for a very 
sub-thermal shock (e = 0.1 so that F ~ 2) with k = 0.05, instability is possible only if Mjso < 6. 
Magnetic fields, which we have not considered, also tend to stabilize the clouds against collapse. 



Gravitational stability of a shocked layer 

We turn now to a more detailed study of the gravitational instability of a shocked cloud, taking 
into account its geometry. The passage of a shock through a cloud that has an approximately 
uniform density (i.e., n > 4) flattens it in the direction of shock propagation. The gravitational 
stability of a pressurized, isothermal sheet of gas was analyzed by Elmegreen Sz. Elmegreen (1978). 
They showed that such a sheet is gravitationally unstable and fragments to form clumps, which may 
be filamentary. When the incident shock is very strong, the sheet is confined by a strong external 
pressure and exhibits effectively incompressible behavior. In this case, the most unstable mode has 
a wavelength that scales with the disk thickness and can be much shorter than the Jeans length 
(Elmegreen & Elmegreen 1978). Such a perturbation can have only a modest growth because the 
total mass is smaller than the critical Bonnor-Ebcrt mass. Only longer wavelengths, which grow 
more slowly, can contain more than a Bonnor-Ebert mass and lead to star formation. Three- 
dimensional numerical calculations by Umekawa (2002) showed that the stable clumps or filaments 
formed by the fastest growing modes can merge with themselves owing to mutual gravitational 
attraction, resulting in formation of more massive spherical cores (see also Umekawa et al. 1999); 
this is presumably equivalent to the growth of longer wavelength instabilities. 

We assume that the mass of the pre-shock cloud is significantly less than the Bonnor-Ebert 
mass, so that it has approximately constant density — i.e., n ^ oo. Correspondingly, its tempera- 
ture, and therefore sound speed, are approximately constant, both before and after the passage of 



-39- 



^=(1 + 737;^) > (50) 



the shock. Elmegreen Sz Elmegreen (1978) analyzed the gravitational instability of an isothermal, 
planar layer of gas. They showed that the parameter that governs the properties of the instability 
is 

(7r/2)GSf. 

where Si is the surface density of the shocked sheet and Pi is the pressure in the shocked intercloud 
medium and therefore at the surface of the shocked cloud. 

The pressure in the midplane of the shocked cloud, Pmid; is greater than that at the surface, 
Pi, due to self gravity: Pi = (1 — ^^)Pmid (Elmegreen Sz Elmegreen 1978). We now show that 
shocked clouds generally have ^ ^ 1, so that we can approximate the shocked sheet as isobaric. 
The surface density of the initial cloud, which is assumed to be uniform, is given in terms of /xq 
(the value of for the pre-shock cloud) by 

So = 0.83 i^-^j nl^'. (51) 

The shock compresses the cloud in the transverse direction by a factor / (which is ~ 0.8 for 
7 = 1.1 and x = 100), so that Si = /"^Sq. Approximating the pressure of the shocked cloud by 
Pi ~ M^Po, we then find that 

Pi /4M2 ^ ^ 

For M > 3 and //q < 1, this is large; as a result, the parameter A is small and the cloud is 
approximately isobaric. 

Elmegreen &; Elmegreen (1978) showed that not all unstable modes lead to gravitational col- 
lapse. For that to occur, two conditions must be satisfied: First, the mode must be large enough 
to capture a Bonnor-Ebert mass, which requires that the wavenumber k satisfies 

kH < 0.91^^/2(1 - ^2)1/4 ^ 0.91^^/2^ (53) 

where H = Cc,i(27rGpinid)'~^/2 is the gravitational scale height of the shocked cloud. The corre- 
sponding e-folding time for the instability is 

tg ~ {A7rGp^,dkHA)-^/\ (54) 

Second, the mode must grow fast enough that it becomes nonlinear during the age of the shock. 
Our results show that the density of the shocked cloud drops by a factor of 2 in 1.5 — 2icc for 
10 < X ^ 100, so we require that the growth time tg be less than 2tcc- In units of the cloud- 
crushing time, 

1/3 
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the growth time of an unstable mode that encompasses at least a Bonnor-Ebert mass is then 

The condition for gravitational collapse to be possible is tg < 2tcc as discussed above. Setting 
F = 1 in equation (43), this leads to 

7/9 

M<7-^, (57) 

so that only shocks of moderate strength can initiate collapse. As discussed above, when the non- 
thermal motions behind the shock are taken into account, the maximum Mach number for collapse 
may be reduced further. The conclusion that shocked clouds are stable against collapse unless the 
Mach number is small differs from that of Vanhala &; Cameron (1998) since they considered a cloud 
that is initially highly centrally concentrated, whereas we have assumed that the initial cloud is 
not self-gravitating and is therefore approximately uniform. 



10. Conclusion 

We have studied the interaction of shock waves with interstellar clouds, using a second-order, 
axisymmetric hydrodynamic code based on the Godunov method with local adaptive mesh refine- 
ment. This problem has been studied previously by KMC, who assumed that there is a discon- 
tinuous density jump between the cloud and the intercloud medium. In reality, one expects that 
the density distribution across the cloud boundary will be smooth due to thermal conduction, non- 
equilibrium cooling effects, or self-gravity (Kornreich &; Scalo 2000). Here we have generalized the 
treatment of KMC to include the effects of smooth cloud boundaries on the shock-cloud interaction. 
We have considered a steady planar shock that strikes an initially spherical cloud with a density 
distribution that is fiat near the center and decreases with radius as a power-law in the envelope. 
Our model can be specified by four numbers: n, the power-law index of the density distribution in 
the envelope; X; the density contrast between the cloud center and the intercloud medium; M, the 
Mach number of the shock; and 7, the ratio of specific heats. 

Our numerical simulations have shown that the clouds are completely destroyed, primarily 
by the Kelvin-Helmholtz instability, irrespective of the initial density distributions, although the 
smoothness of the initial density distribution affects the timescale of cloud destruction. For a cloud 
with n = 2 (the smoothest case we considered), the cloud destruction time is about 3 times longer 
than that of a uniform cloud. The timescale of the cloud destruction is closely related to the 
formation of a density discontinuity (the slip surface) at the upstream side of the shocked cloud. 
The slip surface is subject to the Kelvin-Helmholtz and Rayleigh-Taylor instabilities, and therefore 
significant cloud destruction starts there. Since the growth rates for these instabilities are larger for 
perturbations with shorter wavelengths in the linear regime and since the growth of perturbations 
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with wavelengths shorter than the thickness of the shear layer is greatly suppressed, the growth of 
the instabilities (and accordingly cloud destruction and mixing) are impeded prior to the formation 
of the slip surface. It therefore takes more time for shocks to destroy clouds with smoother density 
gradients. 

We have constructed simple analytic models describing the time evolution of the cloud velocity 
and the circulation that agree well with our numerical results. The vorticity is produced by the 
baroclinic term in the vorticity equation, and can be divided into four components: (1) Tghock; 
which is generated by the initial shock-cloud interacion, (2) Fpost, which is continuously produced 
by the interaction between the shocked cloud and post-shock intercloud flow, (3) Fring, which 
is produced in the intercloud region by the supersonic vortex ring, and (4) Tdoud) which is the 
circulation produced in the cloud and is smaller than other components by a factor of x"^*^^- 
smooth clouds with n = 2, the first component is generated throughout the cloud envelope after 
the initial shock passage, as discussed by Kornreich & Scalo (2000). However, most of the vorticity 
generated inside the cloud shifts towards the upstream side of the cloud, being concentrated in the 
slip surface, and eventually is converted into random motions of small fragments that are detached 
from the main cloud. 

Vortical motions due to hydrodynamic instabilities produce substantial velocity dispersions 

in the shocked cloud. The vortical motions generated in the cloud are converted into random 
motions of small fragments shredded from the parent clouds. Over the wide range of initial model 
parameters considered here, the one-dimensional velocity dispersions are typically 5v r-^ O.lvh, which 
is estimated as about 2 km for T ~ 70 K, % ~ 10^, and M ~ 3. Such velocity dispersions are 
consistent with the observed internal motions of about 1.5 km s~^ in the CNM (Heiles & Troland 
2003) . Recently, Heiles & Troland (2003) suggest that the CNM is in thin sheets with large aspect 
ratios, by assuming that the thickness along the linc-of-sight is estimated to be 6 = iVni/^T-Hl- If we 
apply the same procedure to estimate the aspect ratios, then the shocked clouds torn up into small 
fragments also have large aspect ratios, and therefore, some of the observed CNM clouds could be 
explained by the shocked clouds, instead of thin sheets. 

At the early stage of cloud destruction, the small-scale fiuctuations are dominant due to the 
nonlinear growth of the Kelvin-Helmholtz instability. Thus, the linewidth-size relation in a shocked 
cloud is expected to be relatively fiat. Thereafter, the small-scale fluctuations tend to damp grad- 
ually with time, leading to a linewidth that increases with size. 

Under some circumstances, shock compression may induce gravitational collapse, leading to 
triggered star formation. We have discussed how the critical Bonnor-Ebert mass changes due 
to shock passage. The critical Bonnor-Ebert mass can decrease only when 7 < 4/3. Radiative 
shocks are often assumed to be isothermal, but in fact they are generally sub-isothermal. We 
have determined the density jump across a sub-isothermal shock and shown that, in the absence 
of significant nonthermal motions behind the shock, star formation can be induced in an initially 
uniform cloud by a shock only if it is of moderate strength. Our results indicate that the postshock 
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gas has significant nonthermal motions, and this can further hmit the onset of gravitational collapse. 

We would like to thank M. R. Krumholz for making his valuable IDL codes available to 
us and for stimulating discussion on interstellar turbulence. We also thank R. K. Crockett, D. 
Gies, and J. Walters for a number of valuable conversations. F. N. acknowledges the support 
of the JSPS Postdoctoral Fellowship for Research Abroad. He is also grateful to the Astronomy 
Department, University of California at Berkeley for the kind hospitality during his stay. CFM 
gratefully acknowledges the support of the Miller Institute for Basic Research and of the NSF 
through grant AST-0098365. Part of this work was supported by NASA ATP grant NAG5-12042 
(RIK and CFM) and under the auspices of the US Department of Energy at the Lawrence Livermore 
National laboratory under contract W-7405-Eng-48 (RIK and RTF). Computations were made 
possible by the NSF San Diego Supercomputer Center through NPACI program grant UCB267 
(RIK et al.). 

REFERENCES 

Bacmann, A., Andre, P., Puget, J.-L., Abergel, A. Bontemps, S., Ward-Thompson, D. 2000, A&A, 
361, 555 

Bedogni, R. k Woodward, R. R., 1990, A&A, 231, 481 
Burkert, A. & O'Deh, C. R. 1998, ApJ, 503, 792 
Caselh, P. & Myers, P. C. 1995, ApJ, 446, 665 

Castor, J. I. 1993, in ASP Conf. Ser. 35, Massive Stars: Their Lives in the Interstellar Medium, ed. 
J. P. Cassinelh & E.B. Churchweh (San Francisco: ASP), p. 297 

Chandrasekhar, S. 1961, Hydrodynamic and Hydromagnetic Stability (New York:Dover) 

Ciofh, D. F., McKee, C. F., & Bertschinger, E., 1988, ApJ, 334, 252 

Cox, D. P. & Smith, B. W., 1974, ApJ, 189, L105 

Cowie, L. L. & McKee, C. F., 1977, ApJ, 211, 135 

Dai, W., & Woodward, P.R. 1994, ApJ, 436, 776 

Dame, T. M., Elmegreen, B. G., Cohen, R. S. & Thaddeus, P. 1986, ApJ, 305, 892 
Danforth, C. W., Blair, W. P., & Raymond, J. C. 2001, AJ, 122, 938 
Draine, B. T. & McKee, C. F. 1993, ARA&A, 31, 373 
Elmegreen, B. G., & Elmegreen, D. B. 1978, ApJ, 220, 1051 



-43- 



Elmegreen, B. G., & Lada, C. J., 1977, ApJ, 214, 725 
Elmegrccn, B. G., & Scab, J., 2004, ARA&A, 42, 211 

Falgaronc, E., Lis, D. G., Phillips, T. G., Bouquet, A., Porter, D. H., & Woodward, P. R. 1994, 
ApJ, 436, 728 

Foster, P. N. & Boss, A. P. 1996, ApJ, 468, 784 

Gregori, G., Miniati, P., Ryu, D. & Jones T.W. 2000, ApJ, 543, 775 

Heathcote, S. R. & Brand, P. W. 1983, MNRAS, 203, 67 

Heiles, G. & Troland, T. H. 2003, ApJ, 586, 1067 

Hennebelle, P., Whitworth, A. P., Goodwin, S. P., & Andre, P. MNRAS, 340, 870 
HcnncbcUc, P., Whitworth, A. P., Cha, S.-H., & Goodwin, S. P. MNRAS, 348, 687 
Klein, R. 1. 1999, Jour. Comp. Appl. Math. 109, 123 

Klein, R. I., McKee, G. P., & Golella, P. 1990, in The Evolution of the Interstellar Medium, ed. L. 
Blitz (San Francisco: ASP), p.ll7 

Klein, R. I., McKee, G. F., & Golella, P. 1994, ApJ, 420, 213 

Klein, R. I. &; McKee, G. F. 1994, in Numerical Simulations in Astrophysics: Modeling the Dy- 
namics of the Universe, eds. J. Franco, S. Lizano, L. Aquilar, and E.Daltabuit (Gambridge 
University Press), p. 251 

Klein, R. I. k Wood, D. T. 1998, ApJ, 497, 777 

Klein, R. I., Budil, K. S., Perry, T. S., & Bach, D. R. 2000, ApJS, 127, 379 
Klein, R. I., Budil, K. S., Perry, T. S., & Bach, D. R. 2003, ApJ, 583, 245 
Kornreich, P. k Scab, J. 2000, ApJ, 531, 366 
Koyama, H. k Inutsuka, S. 2002, ApJ, 564, L97 

Landau, L. D., k Lifshitz, E. M. 1959, Fluid Mechanics (Reading: Addison- Wesley) 
Larson, R. B. 1981, MNRAS, 194, 809 
Larson, R. B. 1985, MNRAS, 214, 379 
Larson, R. B. 2005, MNRAS, 359, 211 

Mac Low, M.-M., McKee, C. P., Klein, R. I., Stone, J. M., k Norman, M. L. 1994, ApJ, 433, 757 



-44- 



Mac Low, M.-M. & Klessen, R. 2004, Rev. Mod. Phys., 76, 125 
McKee, C. F. & Cowie, L. L. 1975, ApJ, 195, 715 

McKee, C. F., Hollenbach, D. J., Seab, C. G., & Tielens, A. 1987, ApJ, 318, 674 
McKee, C. F. & Holliman II, J. H. 1999, ApJ, 522, 313 
McKee, C. F. & Ostriker, J. R 1977, ApJ, 218, 148 

Myers, P. C. in Interstellar Processes, ed. D.J. Hollenbach & H. A. Tlironson (D. Reidel Publishing 
Company, Dordrecht), p. 71 

Nittman, J., Falle, S., & Gaskell, P. 1982, MNRAS, 201, 833 

Nulsen, P. E. J. 1982, MNRAS, 198, 1007 

O'Deh, C. R., Henney, W. J., & Burkert, A. 2000, AJ, 119, 2910 

Patnaude, D. J., Fesen, R. A., Raymond, J. C, Levenson, N. A., Graham, J. R., k, Wallace, D. J. 
2002, AJ, 124, 2118 

Poludnenko, A. Y., Prank, A., & Blackman, E.G. 2002, ApJ, 576, 832 

Robey, H. F., Perry, T. S., Klein, R. I., Kane, J. O., Greenough, J. A., h Boehly, T. R. 2002, 
Phys. Rev. Lett., 89, 5001 

Ryutov, D., Drake, R. P., Kane, J., Liang, E., Remington, B. A., & Wood-Vasey, W. M. 1999, 
ApJ, 518, 821 

Schiano, A.V.R., Christiansen, W.A., & Knerr, J.M. 1995, ApJ, 439, 237 
Sgro, A. G. 1975, ApJ, 197, 621 
Spitzer, L. 1982, ApJ, 262, 315 

Stanimirovic, S., k. Heiles, C. 2005, astro-ph/0505106 
Stone, J. M. & Norman, M. L. 1992, ApJ, 390, L17 

Sytine, I. V., Porter, D. H., Woodward, P. R., Hodson, S. W., & Winkler, K.-H. 2000, Jour. Comp. 
Phys., 158, 225 

Tenorio-Taglc, G. & Rozyczka, M. 1986, A&A, 155, 120 

Truelove, J. K., Klein, R. I., McKee, C. F., Holloman II, J. H., HoweU, L. H., Greenough, J. A., & 
Woods, D. T. 1998, ApJ, 495, 821 

Umekawa, M., Matsumoto, R., Miyaji, S., k, Yoshida, T. 1999, PASJ, 51, 625 



-45- 



Umekawa, M. 2002, in The Proceedings of the lAU 8th Asian-Pacific Regional Meeting, Volume II. 
eds. S. Ikeuchi, J. Hearnshaw, and T. Hanawa, (the Astronomical Society of Japan) p. 227 

Vanhala, H. A. T. & Cameron, A. G. W. 1998, ApJ, 508, 291 

Vazquez-Semadeni,E., Ostriker, E. C, Passot, T., Gammie, C. F., & Stone, J. M. 2000, in Protostar 
and Planets IV, eds. V. Mannings, Boss, A. P., k. Russell, S. S. (The University of Arizona 
Press), p. 3 

Woodward, P. R. 1976, ApJ, 207, 484 

Xu, J. & Stone, J. M. 1995, ApJ, 454, 172 



This preprint was prepared with the AAS IM^^X macros v5.2. 



-46- 



Table 1. Model Parameters 



Model 


n 


X 


M 


ao/vco 


7^ 




(Pcl) IPiO 


Resolution 


7 = 5/3, sphere 


ASl 


2 


10 


1.5 


2.61 


6.58 


1.22 


2.91 


-R120 


AS2 


2 


10 


10 


2.61 


6.58 


1.22 


2.91 


jy ID ID TD JD 

-n-SO ) -n-60 > -"-120 ) -"'240 ) -"-480 


ASS 


2 


10 


100 


2.61 


6.58 


1.22 


2.91 


^120 


AS4 


2 


100 


10 


8.50 


34.7 


1.02 


3.58 


-Rl20 


AS5 


4 


10 


10 


1.49 


1.91 


0.611 


4.01 


-^120 


AS6 


4 


100 


10 


2.23 


2.67 


0.510 


8.57 


-^120 


AS7 


8 


10 


1.5 


1.15 


1.23 


0.306 


5.62 


Rl20 


ASS 


8 


10 


10 


1.15 


1.23 


0.306 


5.62 


R30 , -^60 , -^120 , -R240 , -R48O , -^960 


AS9 


8 


10 


100 


1.15 


1.23 


0.306 


5.62 


-R12O 


ASIO 


8 


10 


1000 


1.15 


1.23 


0.306 


5.62 


R12O 


ASH 


8 


100 


10 


1.27 


1.28 


0.255 


23.0 


-R12O 


AS12 


16 


10 


10 


1.04 


1.07 


0.153 


7.22 


R12O 


AS13 


24 


10 


10 


1.02 


1.03 


0.102 


7.97 


R12O 


ASM 


oo 


10 


10 


1.00 


1.00 


0.00 


10.0 


R30 1 R-GO ) R12O , -^240 


AS15 


oo 


10 


100 


1.00 


1.00 


0.00 


10.0 


^120 


AS16 


oo 


100 


10 


1.00 


1.00 


0.00 


100.0 


R12O 


7 = 1.1, 


sphere 














ISl 


2 


10 


10 


2.61 


6.58 


1.22 


2.91 


Rl2Q 


IS2 


8 


10 


10 


1.15 


1.23 


0.306 


5.62 


R12O 


7 = 5/3, cylinder 


ACl 


2 


10 


10 


2.75 


2.78 


1.22 


3.47 


R12O 


AC2 


8 


10 


10 


1.25 


1.11 


0.306 


6.59 


R12O 



Note. — In the eighth column, (pci) denotes the mean cloud velocity. In the sixth column, TZ 
denotes the cloud mass normalized with mo = 4TTr^Pco/S and Trr^oPcO for spherical and cylindrical 
clouds, respectively. To estimate the cloud mass of a cylindrical cloud, we set the vertical height to 
be equal to the core radius Tco- In the last column, i?„ means a grid resolution of n cells per core 
radius measured at the finest level of refinement (level 3). 
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Table 2. Dependence on n (x = 10, M = 10, R^o) 



iViOQGl 




AS2 




AS5 




A 


AS12 


AS13 


A m /I 


P cir 3,111 Gt G r 


n = 2 


n = 4 




r) — R 
fl — o 


n = 16 


n = 24 




'-slip/^cc 




4.91 




1.98 




1.71 


0.722 


0.335 


0.00 


^drag/^cc 




7.51 




3.75 




3.26 




3.15 




2.90 


2.52 


^dest / ^cc 




10.2 




7.40 




5.91 




5.34 




4.47 


3.45 


^mix / ice 




4.73 




8.04 




6.59 




6.31 




6.17 


4.91 


tm /ice 




12.9 




5.96 




4.82 




4.54 




4.24 


3.80 


a/rco 


3.04 


(3.08) 


2.07 


(1.98) 


1.73 


(1.66) 


1.75 


(1.75) 


1.76 


(1.75) 


2.21 




3.45 


(2.67) 


1.85 


(1.05) 


1.48 


(1.29) 


1.62 


(1.51) 


1.79 


(1.75) 


1.63 


c/a 


1.12 


(0.865) 


0.893 


(0.529) 


0.855 


(0.775) 


0.928 


(0.862) 


1.02 


(0.998) 


0.738 


{P)/{P0) 


2.55 


(1.59) 


2.20 


(1.97) 


1.90 


(1.84) 


1.65 


(1.59) 


1.48 


(1.42) 


1.24 


(Vz) /Vh 


0.173 


(0.124) 


0.182 


(0.186) 


0.200 


(0.206) 


0.173 


(0.170) 


0.162 


(0.162) 


0.130 




0.0758 


(0.0840) 


0.0953 


(0.0788) 


0.102 


(0.0840) 


0.113 


(0.105) 


0.116 


(0.113) 


0.115 


SVz jv-b 


0.0990 


(0.0703) 


0.0870 


(0.0701) 


0.158 


(0.119) 


0.164 


(0.159) 


0.181 


(0.176) 


0.182 


-T/rcoVb 


2.38 




0.501 




1.25 




1.20 




1.17 




1.50 



Note. — All the timc-dcpcndcnt quantities arc estimated at t = tm- The initial cloud mass md.o is adopted as the threshold 
mass in computing the global quantities, except that the values in parentheses are measured with respect to the initial core 
mass, rricofi- Note that the radii a and c are normalized by rco- If the radius a is normalized by its initial value ao, then the 
radius of the shocked cloud decreases monotonically with n. See §4.3. 



Table 3. Dependence on n (x = 100, M = 10, i?i2o) 



Model 




AS4 




AS6 




ASH 


AS16 


Parameter 




n = 2* 


n = 4 




n = 8 


n = GO 


^slip/^cc 




3.50 




1.63 




1.10 


0.00 


^drag /^cc 




10.5 




5.22 




5.01 


3.59 


^dest/^cc 




12.1 




6.51 




6.50 


3.02 


^mix/^cc 




12.0 




6.24 




5.69 


4.72 


tm /ice 




> 15.0 




7.05 




6.57 


4.01 


a/Vco 


8.14 


(0.385) 


3.42 


(1.91) 


3.16 


(2.83) 


3.15 


C/Vco 


26.6 


(4.45) 


15.3 


(3.41) 


11.2 


(9.50) 


3.15 


c/a 


3.27 


(11.5) 


4.48 


(1.79) 


3.57 


(3.35) 


2.65 


(P> / (P0> 


2.01 


(0.240) 


0.959 


(0.390) 


0.245 


(0.289) 


2.65 




0.239 


(0.286) 


0.181 


(0.314) 


0.184 


(0.217) 


0.240 




0.0571 (0.0582) 


0.0810 


(0.0418) 


0.102 


(0.0840) 


0.0987 


SVz /Vb 


0.142 


(0.0989) 


0.0870 


(0.0701) 


0.175 


(0.106) 


0.236 


—T/rcoVb 


17.0 




10.4 




8.03 




6.51 



The time-dependent quantities are evaluated at t = 16.2 tec. 



Note. — The time-dependent values are estimated aX t = tm- The initial cloud mass mci,o is adopted as the threshold mass 
to compute the global quantities, except in the values in parentheses which are the global quantities measured in the initial 
core mass, mco, • 
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Table 4. Dependence on M (x = 10, -R120) 



Parameter n = 8^ n = 2 



Model 




AS7 




AS9 




ASIO 




ASl 




AS2 




AS3 


M 




1.5 




100 




1000 




1.5 




10 




100 






0.417 




0.750 




0.750 




0.417 




0.742 




0.750 


^drag/^cc 




17.2 




3.25 




3.25 




21.4 




7.51 




6.92 


^dest/^cc 




12.2 




5.88 




5.87 




32.8 




10.2 




10.1 


^mix/^cc 




15.8 




6.52 




6.60 




25.1 




4.73 




4.92 


im /ice 




7.51 




5.15 




5.13 




> 32.7 




12.9 




13.6 


a/r-co 


1.82 


(1.74) 


1.80 


(1.71) 


1.80 


(1.78) 


3.89 


(4.02) 


3.04 


(3.08) 


3.18 


(3.24) 




1.49 


(1.18) 


1.60 


(1.48) 


1.60 


(1.50) 


5.20 


(4.39) 


3.45 


(2.67) 


3.46 


(2.32) 


c/a 


0.816 


(0.680) 


0.889 


(0.864) 


0.888 


(0.845) 


1.34 


(1.09) 


1.12 


(0.865) 


1.09 


(0.716) 


{P)/{P0) 


1.33 


(1.37) 


1.86 


(1.70) 


1.86 


(1.72) 


1.25 


(0.712) 


2.55 


(1.59) 


2.55 


(1.50) 


(Vz) /Vh 


0.193 


(0.197) 


0.182 


(0.193) 


0.182 


(0.193) 


0.114 


(0.0904) 


0.173 


(0.124) 


0.159 


(0.125) 




0.0629 


(0.0508) 


0.104 


(0.0900) 


0.104 


(0.0908) 


0.0475 (0.0410) 


0.0758 


(0.0840) 


0.0781 


(0.0732) 




0.0879 


(0.0646) 


0.172 


(0.146) 


0.164 


(0.140) 


0.0738 (0.0654) 


0.0990 


(0.0703) 


0.101 


(0.0711) 


-T/rcoVb 


1.55 




1.22 




1.22 




2.23 




2.38 




2.37 





^see Table 2 for the global quantities of x = 10 

Note. — The time-dcpcndcnt values arc estimated at t = tm- The initial cloud mass m^.o is adopted as the threshold mass 
to compute the global quantities, except in the values in parentheses which are the global quantities measured in the initial 
core mass, mco, ■ 



Table 5. Radiative Clouds (x = 10, M = 10, i?i2o) 



Model 




ISl 




IS2 


Parameter 


n = 2 


n = 8 


^drag/'-cc 




0.500 




1.56 


^dcst / ^cc 




7.23 




6.71 


^niix / ^rc 




6.10 




5.56 


i„./''c-i- 




8..S7 




1.9<) 


a/r-co 


2.41 


(1.29) 


1.11 


(0.830) 


c/r-co 


6.74 


(8.19) 


4.54 


(3.20) 


c/a 


2.79 


(6.36) 


1.15 


(3.86) 


ip) 1 (Po) 


10.8 


(6.12) 


7.05 


(5.65) 


{Vz) /Vb 


0.129 


(0.0649) 


0.164 


(0.196) 


Sv.a,/Vb\ 


0.0541 


(0.0500) 


0.0553 


(0.0532) 


SVz /vt, 


0.0989 


(0.149) 


0.106 


(0.0793) 


-T/rcoVb 


0.140 




0.327 





Note. — We note that the circulation associated with the supersonic vortex ring, normalized with rcoffc, contributes more for 
radiative clouds than for nonradiative clouds. Since this has the opposite sign from the other contributions to the circulation, 
the values for the total circulation listed above are smaller than those in Table 2. However, the peak values of the circulation 
are close to those of the nonradiative cases. 
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Table 6. Cylindrical Clouds (x = 10, M = 10, R^o) 



Model 




ACl 




AC2 


Parameter 


n = 2=- 


n = 8 


^drag/^cc 




2.04 




2.38 






17.1 




6.25 


i mix /ice 




18.7 




11.3 


im /tec 


> 19.0 




9.14 


a/fco 


3.43 


(3.53) 


2.23 


(2.23) 


c/rco 


3.00 


(2.84) 


3.24 


(3.15) 


c/a 


0.876 


(0.805) 


1.45 


(1.41) 


(P) / (Po) 


2.46 


(2.00) 


1.50 


(1.37) 


{v.) /vt 


0.0915 


(0.0975) 


0.111 


(0.111) 


Svza/Vb 


0.0758 


(0.0635) 


0.0893 


(0.0893) 


SVz /vi, 


0.0885 


(0.0662) 


0.136 


(0.136) 


-T/rcoVb 


2.03 




2.24 





^For n = 2, the time-dependent global quantities are evaluated at t = 12.9 tec, which is the maximum expansion time for the 
spherical case. 



Note. — The time-dependent values are estimated 3.t t = tm- The initial cloud mass mci,o is adopted as the threshold mass 
to compute the global quantities, except in the values in parentheses which are the global quantities measured in the initial 
core mass, mco, . 
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Fig. 1. — Results of convergence studies for (a) {n,x,M) = (8,10,10) and (b) (2, 10, 10). The 
ordinate denotes the relative error, while the abscissa denotes the number of grid cells per radius, 
which is a measure of the spatial resolution. See the text for the definition of the relative errors. 
The global quantities are compared at t = 3 tec and 6 tec for n = 8 and 2, respectively. 
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Fig. 2. — Snapshots of the density distribution for model ASS with n = 8, x = 10, and 7 = 5/3 
interacting with a strong shock of M = 10 at seven different times: (a) initial state {t = —0.39 tec); 

(b) 1.24 tec, (c) 2.00 tec, (d) 4.00 tec, (e) 6.00 tec, (f) 8.00 tee, and (g) 10.00 tec, where the evolution 
time is measured from the stage at which the shock has reached at z/rco = — 1 on the 2;-axis. Tha 
abscissa and ordinate denote the ru- and z-coordinates, respectively, normalized by the core radius, 
Tco- The initial planar shock was initiated at z/rco = —2.1. The arch-like density jump appearing 
at z/rco ~ —1 and —0.85 in panels (b) and (c), respectively, is a bow shock, which is formed by the 
shock reflected at the cloud surface. The black spot appearing at (cc^/rco, z/rco) ~ (0-2,5) in panel 

(c) is a supersonic vortex ring, which propagates, with almost the same speed as the intercloud 
shock, up to {'Uj/reo,z/reo) — (0.5,11.4) and (0.6,17.6) in panels (d) and (e), respectively. 
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Fig. 3. — Snapshots of the density distribution for model AS2 with n = 2, x = 10, and 7 = 5/3 
interacting with a shock of M = 10 at seven different times: (a) initial state (t = —5.76 ^cc), (b) 
1.93 tec, (c) 3 tec, (d) 6 tec, (e) 9 tec, (f) 12 tec, and (g) 15 tee, where the evolution time is measured 
from the stage at which the shock has reached at z/vco = —1 on the z-axis. The initial planar 
shock was initiated at z/vco = —18. 
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Fig. 4. — (a), (b), and (c): Snapshots of the density distribution for model AS4 with {n,x,M) = 
(2,100,10) at (a) t = 3tcc, (b) t = 6tcc, and (c) t = 9tcc- (d), (e), and (f): Snapshots of the 
density distribution for model ASH with {n,x,M) = (8, 100, 10) at (d) t = 2tcc, (e) t = 4tcc, and 
(f) i = 6tcc- 



4:.? 4:j C.I ci 




Fig. 5. — Mach scaling. Comparison of the density distributions between models ASS and AS9 
with n = 8 and X = 10: (a) model AS9 (M = 100) at t = 3.5 tec, (b) model ASIO (M = 1000) at 
t = 3.5tcc , (c) model AS9 (M = 100) at t = 6.4tcc, (d) model ASIO (M = 1000) at t = 6.4tcc- 

The morphology of the cloud is remarkably similar at the first time, whereas small differences in 
the small-scale structures begin to appear at the second time. Note that the global quantities agree 
very well even at the later times (see the text for more detail). 
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Fig. 6. — (a) and (b): Time evolution of the cloud shape and the cloud velocity for two different 
models with {n,x,M) = (8,10,10) {thick curves) and (oo,100, 10) {thin curves). Panel (a) shows 
the time evolution of the rms radial radius a {solid curves) and axial radius c {dotted curves). 
Panel (b) shows the time evolution of the mean cloud velocity {solid curves). For comparison, the 
analytic solution is indicated by dotted curves, (c) and (d):S&me as the upper panels but for the 
models with (n, x,-^^) = (2, 10, 10) and (4, 10, 10). (e) and (f):Swa.e as the upper panels but for 
the models with 
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Fig. 7. - Distributions of (a) density, (b) Vr, and (c) Vz, along r = 0.5rco, for n = 8, x = 10, and 
M = 10. Solid, dotted, dashed, and dash-dot curves indicate the profiles at t = 0.249tcc) 1-31 tec, 
2.13 icc) and 3.04 respectively. Distributions of (d) density, (e) Vr-, and (f) Vz, along r = O.Srco, 
for n = 2, X = 10, and M = 10. Sohd, dotted, dashed, dash-dot, and dash-dot-dot-dot curves 
indicate the profiles at t = 1.49 tco 2.82 tec, 4.05 tec, 5.18 tco and 6.41 tco respectively. 
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Fig. 8. — Time evolution of the radial (dotted curve) and axial (solid curve) velocity dispersions 
for (a) n = 8, X = 10, and M = 10 and (b) n = 2, x = 10, and M = 10. 




Fig. 9. — Time evolution of fractional masses with densities greater than 2 pii, normalized by the 
initial cloud mass. 
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Fig. 10. — Evolution of the mass fraction for the models with (n,x,M) = {8, 10, 10) at times (a) 
t = 0, (b) 1.06 tec, (c) 6.04 tec, and (d) 8.00 tec, and for the models with (n,x,M) = (2, 10, 10) at 
times (e) t = 0, (f) 1.40 tco (g) 7.28 ico and (h) 29.5 tec- The histograms indicate the fraction of 
mass contained within a corresponding density bin with a width of 0.2 pco, which is normalized by 
the total cloud mass mci,o- 




Fig. 11. — Time evolution of circulations for clouds with (a) x = 10 and clouds with (b) x = 100. 
Thick solid, thin solid, and dashed curves denote the time evolution of the total circulation for the 
models with n = 2, 8, and oo, respectively. 
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Fig. 12. — Density and vorticity distributions for the model with (n, Xi = (8, 10, 10). Panels (a) 
and (b): the density distribution and vorticity distribution at t = 0.930 tec- Panels (c), (d), and 
(e): the density distribution, vorticity distribution, and the velocity distribution at t = 2.86 tec- 
The velocity vectors are measured in the frame that moves with the postshock flow. Panel (f): 
cumulative circulation at t = 2.86 ice- 



Fig. 13. — Density and vorticity distributions for the model with (n,x, M) = (2,10,10). Panels 
(a) and (b): the density distribution and vorticity distribution at t = 2.02tcc- Panels (c), (d), and 
(e): the density distribution, vorticity distribution, and the velocity distribution at t = 6.65 tec- 
The velocity vectors are measured in the frame that moves with the postshock flow. Panel (f): 
cumulative circulation at t = 6.65 tec- 




Fig. 14. — Evolution of the mass ratio for the models with (n, x, M) = (8, 100, 10) and (2, 100, 10), 
which are plotted with solid and dotted curves, respectively. 
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Fig. 15. — Synthesized spectra for the model with (n, %, M) = (8, 10, 10) at (a) t = 4.00 tec, (b) 
S.OOtcc, and (c) 13.2tcc, and for the model with (n,x,M) = (2,10,10) at (d) t = 6.00 ice, (e) 
12.0tcc) and (f) 20.1 The spectra are normalized by the maximum at (a) and (d) for n = 8 and 
2, respectively. The histograms denote the mass contained within a corresponding velocity bin with 
a width of O.OlSfji- The abscissa denotes the velocity normalized by the velocity of the postshock 
ambient medium vu. The integration is done over the entire cloud under the assumption that the 
line-of-sight is parallel to the z-axis. 
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Fig. 16. — linewidth-size relation in shocked clouds with (n, x, M) = (8, 10, 10) at three different 
times: (a) t = 6 tec, (b) 8 tec, and (c) 10 ice- This model is followed with the three-dimensional 
simulation with an intermediate resolution i?60- The abscissa denotes the beam size, while the 
ordinate is the velocity dispersion estimated in the beam with a size R. The velocity dispersions 
arc measured in units of the sound speed at the center of the preshock cloud, Cco- Solid curves 
show the velocity dispersions measured from the direction whose line-of-sight is identical to the 
z-axis. Dotted and dashed curves show the velocity dispersions measured from the direction whose 
line-of-sight is perpendicular to the z-axis. For dotted and dashed curves, the beam center is put 
on the center of gravity, zq, and zq + ^co, respectively. 



